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PREFACE 


This  report  summarizes  the  procedure  for  calibrating  thermistors  as  the 
temperature  sensing  elements  in  simulated  skin  sensors.  The  sensors  will  be 
used  in  the  new  Advanced  Thermal  Response  Data  Acquisition  and  Analysis  System 
(ATRES/DAAS) .  This  work  is  being  done  under  the  project  "Thermal  Protection, 
1L162786AH98CABOO" ,  administered  by  the  Physics  &  Engineering  Branch  (PEBr), 
Fiber  and  Polymer  Science  Division  (FPSD),  Soldier  Science  Directorate,  U.S. 
Army  Natick  Research,  Development,  and  Engineering  Center. 

All  experiments  are  conducted  at  the  Physics  &  Engineering  Branch  at  the 
U.S.  Army  Natick  Research,  Development,  and  Engineering  Center  (Natick).  The 
citation  of  trade  naunes  in  this  report  does  not  constitute  official 
endorsement  or  approval  of  use  of  an  item. 
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THERMISTOR  CALIBRATION  PROCEDURE  FOR  SIMULATED  SKIN  SENSORS 


INTRODUCTION 

A  new^  state-of-the-art  data  acquisition/ instrumented  manikin  system  is 
being  developed  by  the  Physics  and  Engineering  Branch  at  Natick  to  study  the 
degree  of  soldier  protection  afforded  by  various  uniforms  and  materials 
against  thermal  insults  such  as  flames,  CO2  lasers,  and  thermonuclear 
weapons^ The  Advanced  Thermal  Response  Data  Acquisition  and  Analysis  System 

«  (ATRES/DAAS)  consists  of  four  units:  1)  remote  data  collection  sites,  2) 

master  data  collection  control,  conversion,  and  formatting,  3)  data  processing 

•  and  storage,  and  4)  graphical  analysis  and  presentation.  This  report  deals 
with  the  remote  data  collection  sites,  more  specifically,  the  temperature 
sensing  elements  that  are  mounted  in  system  manikins. 

Previous  data  acquisition/ instrumented  manikin  systems  used 
thermocouples  as  the  temperature- sensing  elements.  The  new  ATRES/DAAS  system 
uses  thermistors.  The  previous  thermocouple-mounted  sensors  were  problematic 
in  the  low  signal-to-noise  ratio,  which  is  characteristic  of  thermal 
electromotive  forces  (emfs).  The  high  voltage  output  (volts)  that  can  be 
realized  with  thermistors  virtually  eliminates  this  problem. 

The  thermistors,  as  with  the  previous  thermocouples,  are  mounted  inside 
simulated  skin  sensors.  These  sensors  are  disks,  2.5  cm  in  diameter  and  1.25 
cm  thick,  made  from  a  60/40  mixture  (by  weight)  of  alpha-cellulose-filled  urea 
formaldehyde  and  silica.  The  disks  have  thermal  and  optical  properties 
similar  to  those  of  human  skin^.  Skin  burn  severity  can  be  predicted  from  the 
temperature  recorded  over  time  by  the  ATRES/DAAS  system. 

►  This  report  describes  the  development,  implementation,  and  results  of 

the  thermistor  calibration  procedure  currently  being  employed.  First  an 

*  uncertainty  analysis  is  performed  in  order  to  determine  appropriate  parametric 
measurement  procedures  needed  to  obtain  the  required  accuracy  (±  0.5®C).  The 
calibration  procedure,  based  on  the  analysis,  is  then  described  in  detail. 
Finally,  results  to  date  are  presented,  which  show  temperature  measurement 
precision  within  the  predicted  uncertainty  and  required  accuracy. 
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THEORETICAL  BACKGROUND 


Resistance  thermometry  is  a  temperature  measurement  technique  that  uses 
the  bhahge  in  electrical  resistance  of  a  material  as  its  temperature  changes, 
the  two  basic  classes  of  resistance  thermometers  are  resistance  temperature 
detectors  (RTD)  and  thermistors:  the  former  employs  a  conductor  and  the 
latter,  a  semiconductor. 

I 

The  Advanced  Thermal  Response  Data  Acquisition  and  Analysis  System  being 
developed  at  Natick  uses  negative  temperature  coefficient  (NTC)  thermistors  as 
the  temperature-sensing  elements  of  the  simulated  skin  sensors. 


Thermistor  Type 

Thermistors  were  purchased  from  Thermometries^  as  B05PA103N.  The 

numbers  and  letters  have  the  following  significance: 

B  -  Bead  style:  denotes  a  glass-coated  bead. 

05  -  Bead  diameter:  0.005  inch  nominal  diameter. 

P  -  Lead  Configuration:  denotes  opposite  leads. 

A  -  Material  system:  determines  the  optimal  range  of  resistance  sensitivity 
for  a  particular  temperature  range. 

103  -  Zero  power  resistance:  resistance  at  25®C,  10,000  ohms. 

N  -  Tolerance:  denotes  +-  25%  resistance  tolerance  at  25®C. 

BOS  Thermal  and  Electrical  Properties 

Thermal  time  constant:  in  still  air  -  0.12  s 

water  plunge  -  0.005  s 

Dissipation  constant:  in  still  air  -  0.045mW/®C 

water  plunge  -  0.23mW/®C 
Resistance  range:  IK  to  lOM  ohms 

Maximum  power  rating:  0.006  watts. 


Thermistors  have  two  advantages  over  thermocouples.  They  have  a  higher 
signal-to-noise  ratio  (-10^),  and,  because  of  their  smaller  size,  they  respond 
much  faster  to  temperature  change.  However,  thermistors  are  fragile,  have 
highly  nonlinear  resistance/ temperature  characteristics,  are  not 
interchangeable,  and  are  expensive  when  close  tolerance  is  required. 
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Linear  Approximation 


M 


% 


Although  the  resistance/temperature  characteristics  of  thermistors  are 
highly  nonlinear,  a  linear  approximation  is  obtained  when  the  natural  log  of 
thermistor  resistance  is  plotted  versus  the  reciprocal  of  absolute 
temperature.  In  this  case  we  have  the  relation 

lni?j,=A+-4 

T  J. 

where  T  is  the  absolute  temperature  (Kelvin)  and  p  is  a  par2uneter  dependent 
upon  the  material.  Here,  the  slope  P  is  considered  constant  although  this 
approximation  holds  only  over  a  relatively  small  temperature  interval.  For 
temperatures  between  0®C  and  50®C,  interpolation  errors  introduced  by  assuming 
a  constant  P  are  on  the  order  of  0.01®C  over  a  10®C  temperature  span,  0.04®C 
over  a  20®C  span,  0.1*C  over  a  span  of  aO^C,  and  over  a  SO^C  span. 

The  thermistor  resistance  at  a  reference  temperature  is  given  by  the 
expression® 

(2) 

O 


where  is  the  resistance  at  reference  temperature  T^.  By  subtracting  Eg.  2 
from  Eq.  1  the  standard  thermistor  equation  is  obtained® 


(3) 


The  value  of  P  is  dependent  upon  material,  temperature,  and  thermistor 
construction  and  is  considered  unique  for  each  thermistor.  With  two 
temperature/resistance  data  points,  p  can  be  directly  determined  using  Eq.  3. 
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Polynomial  Approximation 


A  inor©  accurat©  approach  us©s  a  standard  curve-fitting  technigue  with 
th©  assumption  that  the  In  Rj  is  a  polynomial  in  l/T^  or  vice  versa 


lni?7=A.+-^  +  -4  +  .  . 

T  O  fjt  J,2 


<4) 


These  expressions  give  accurate  representation  of  temperature/resistance 
characteristics  over  wide  temperature  spans.  Thermometries,  Inc.  reports  an 
excellent  curve  fit  using  a  third-degree  polynomial  with  the  squared  term 
eliminated®.  Equations  4  and  5  become 


lni?i.=Ao+ 


T 


(6) 


and 

1  =3^+31  (lni?i.) +33 

Eliminating  the  squared  term  reduces  the  number  of  points  needed  to  determine 
polynomial  constants  but  introduces  some  error.  For  example,  for  the 
temperature  span  0<T<150‘’C  the  maximum  error  is  0.045‘’C. 

As  previously  stated,  P  is  the  slope  of  the  In  versus  1/T  (absolute) 

curve 

dilnRj.)  (8) 

d(l/T) 

Using  Eg.  6,  the  expression  for  P  is 

(9) 
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Linear  Voltage  Divider 


Thermistor  resistance/temperature  measurements  are  generally  obtained 
through  voltage  measurements.  The  voltage  divider  is  the  simplest  circuit 
used  in  thermistor  networks  (Fig.  1).  In  this  case  the  output  voltage  ej,(T) 
is  measured  across  the  fixed  resistor  R 


Solving  for  Rj  we  obtain 


e<,(D  =e 


(10) 


(11) 


Advantages  to  using  this  circuit  include  increasing  output  voltage  with 
increasing  temperature  and  the  ability  to  include  the  loading  of  any  external 
measurement  circuitry  in  the  value  of  R.  The  value  of  R  is  chosen  such  that 
there  is  maximum  linearity  over  the  temperature  range  of  interest^. 


Fig  1.  Voltage  divider  circuit  for  thermistor  resistance  measurements. 

Thermal  Time  Constant  Considerations 

The  thermal  time  constant  of  a  thermistor  in  air  is  on  the  order  of  0.12 
seconds.  This,  combined  with  the  extreme  thermistor  fragility,  makes  it 
extremely  difficult  to  calibrate  a  thermistor  as-is.  That  is,  in  order  to 
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obtain  the  thermistor  must  be  held  at  the  constant  reference  temperature 

long  enough  for  a  resistance  measurement  to  be  obtained.  With  a  0,12  s 
thermal  time  constant,  the  thermistor  reacts  quickly  to  subtle  temperature 
Changes  due  to  air  currents,  etc.  Therefore,  the  thermistors  are  first 
mounted  inside  the  simulated  skin  sensor.  This  insulates  the  thermistors  from 
slight  temperature  changes  and  protects  them  from  breakage  due  to  handling. 
With  the  thermistors  mounted  inside  the  sensors,  the  time  constant  of 
interest  is  that  of  the  sensor  material  covering  the  thermistor.  The  time  t 
required  for  heat  to  travel  a  distance  X  in  a  material  with  thermal 
diffusivity  K  can  be  approximated  by  the  expression 

t=i£!.  (12) 

K 

The  thermal  diffusivity  K  is  defined  as 

K=^  (13) 

pCp 

where  k  is  the  thermal  conductivity,  p  is  the  material  density  and  Cp  is  the 
specific  heat  at  constant  pressure.  For  a  60/40  mixture  of  urea- formaldehyde 
and  silica,  the  values  of  k  and  pCp  are  reported  to  be® 
k  s  0.00131  cal/cm-s-®C, 
pCp  =  0.65  cal/cm^-®C. 

These  values  give  a  thermal  diffusivity  value  of  K  —  0.0020  cm  /s. 

The  sensors  are  manufactured  with  an  approximately  0.05  cm  thick  of 
urea-formaldehyde  and  silica.  The  time  constant  calculated  for  this  material 
thickness  using  Eg.  12  is  about  1,25  seconds.  The  time  for  the  applied  heat 
to  raise  the  temperature  of  the  entire  disk  is  also  of  interest.  This  will  be 
important  in  calibrating  the  thermistors  because  it  is  important  to  have  some 
idea  how  long  it  will  take  for  a  sensor  to  reach  thermal  equilibrium.  For 
x=0.65  cm,  which  is  about  half  the  thickness  of  the  sensor,  the  thermal  time 
constant  is  computed  to  be  about  5  min. 

This  relatively  long  thermal  time  constant  indicates  the  importance  of 
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allowing  the  skin  sensor  assembly  to  reach  an  equilibrium  temperature  before 
attempting  to  take  temperature/resistance  readings. 


UNCERTAINTY  ANALYSIS 

In  developing  a  thermistor  calibration  procedure  it  is  advantageous  to 
perform  an  uncertainty  analysis  in  order  to  determine  the  most  accurate 
measurement  methods  by  which  to  quantify  the  parameters  involved.  Here  the 
analysis  is  based  on  the  standard  root- sum- squares  method. 

If  it  is  assumed  that  the  uncertainties  (W^,  W2  f.../Wu)  in  the 
independent  variables  (X,,  Xg, are  all  given  with  the  same  probability, 
then  the  uncertainty  (Wp)  in  the  dependent  variable  (R)  is  given  by  the 
expression 


We  are  interested  in  the  accuracy  with  which  we  can  measure  both  absolute  and 
relative  (AT)  temperature. 


Absolute-Temperature  Uncertainty 


The  uncertainty  analysis  is  based  on  the  linear  approximation.  Using 
Eq.  3  and  solving  for  T,  expressions  can  be  obtained  for  absolute-temperature 
as  a  function  of  th©rrnistor  r©sistanc©s 


r= 


r^ln(V^r)+P 


(15) 


For  th©  above  expression,  there  is  expected  uncertainty  in 
By  inserting  Eg*  15  into  Eg.  14  a  probable  estimate  of  the 
obtained; 


11/2 


Pf  Ryf  ^TO* 

uncertainty  in  T  is 

(16) 
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where 


dT  .  -pTo 

dRj.  i?j.  ( Tp  In  {Rj./Rt)  +  P  )  *  ' 


dT  _ _ 

ai?r„  i?r^(r„ln(VV+P)'' 


3r  _ _ gf _ 

ar^  (r5ln(i?j./i?j.^)+p)*' 


(17) 

(18) 

(19) 


^'p  Tq  ^ti{Ry/RxJ^ 

(r^ln(i?i./i?r.)+p)*  ■ 

The  uncertainty  in  is  obtained  by  applying  Eq.  11  to  Eg.  14 
with  expected  uncertainties  in  R,  e^,  and  e^(T): 


where 


J  ^ 


(20) 


(21) 


(22) 


aRr_  R 
a®,  ejT)' 


(23) 


3Rj.  _  Re^  ^24j 

ae<,(r)  (e„(r))*' 

Using  Eqs.  16-24  the  uncertainty  in  Rj  and  the  resulting  uncertainty  in 
T  was  calculated  for  R^^^  tolerances  of  25%,  10%,  5%,  and  1%,  while  varying  the 
tolerance  on  R  for  values  of  5%,  1%,  and  0«1%.  For  this  study  the  P  value  is 
obtained  from  the  manufacturer^  as  well  as  the  expected  deviation^°.  It  is 
assumed  that  the  source— voltage  is  set  using  a  digital  voltmeter  and  the 
output  voltage  is  recorded  using  LABTECH*  NOTEBOOK.  Values  and  tolerances  for 
P,  ©g,  e^(T),  and  T^,  are  given  as 


t 
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P  =  3395  k’, 

W-  =  101.85  k’”, 
e  =  2.5  Volts, 

W®  =  0.01  Volts  (readability  of  digital  voltmeter), 
e®^T)  =  1.32  Volts  @  44®C  (calculated  using  Eqs.  3  and  10), 
wL,t^  =  1.32x10'^  Volts”, 

T®  ^>2500, 

=  0  (tolerance  unknown,  assumed  negligible). 

The  resulting  values  are  given  in  Table  1.  Also  included  in  Table  1  are  the 
uncertainties  associated  with  the  individual  terms.  These  are  helpful  in 
determining  the  major  contributors  to  the  overall  uncertainty  of  R-j-  and  T. 
These  terms  are  designated  as  U  and  are  determined  by  multiplying  the 
parameter  sensitivity  by  its  respective  tolerance  as  follows: 


(25) 


e-(T) 


(26) 


Re^  67 

[eAT)P 


(27) 


rr  P  To 

Rr[r<,ln(i?r/i?r) +PP 


(28) 


p  tI 

[r„ln(VV^P^' 


(29) 


Tq  ln{Rf/R]~J)  ^ 

[r„ln(V«r)+P]^  ^ 


(30) 
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Table  1.  Uncertainty  in  and  absolute  T  with  variation 


in  and  R  tolerances. 


Tolerances 

Uncertainty 

in  R, 

Uncertainty  in  T 

(abs) 

Bio 

R 

% 

Hes 

Seol 

%r 

% 

%To 

h 

25% 

5% 

1394 

152 

3.8 

1402 

8.2 

7.4 

0.6 

11 

10% 

5% 

1394 

152 

3.8 

1402 

8.2 

3.0 

0.6 

8.7 

5% 

5% 

1394 

152 

3.8 

1402 

8.2 

1.5 

0.6 

8.4 

1% 

5% 

1394 

152 

3.8 

1402 

8.2 

0.3 

0.6 

8.2 

25% 

1% 

279 

152 

3.8 

317 

1.8 

7.4 

0.6 

7.6 

10% 

1% 

279 

152 

3.8 

317 

1.8 

3.0 

0.6 

3.5 

5% 

1% 

279 

152 

3.8 

317 

1.8 

1.5 

0.6 

2.4 

1% 

1% 

279 

152 

3.8 

317 

1.8 

0.3 

0.6 

2.0 

25% 

0.1% 

28 

152 

3.8 

155 

0.9 

7.4 

0.6 

7.5 

10% 

0.1% 

28 

152 

3.8 

155 

0.9 

3.0 

0.6 

3.1 

5% 

0.1% 

28 

152 

3.8 

155 

0.9 

1.5 

0.6 

1.8 

1% 

0.1% 

28 

152 

3.8 

155 

0.9 

0.3 

0.6 

1.1 

A 

B 

C 

D 

E 

F 

6 

H 

I 

J 

As  would  be  expected/  the  overall  uncertainty  in  column  F)  and 

T  (Wy,  column  J)  decreases  as  the  and  R  tolerances  become  smaller.  From 
this  table  it  is  evident  that  the  R  tolerance  is  critical  with  as  much  as  an 
86%  decrease  in  temperature  uncertainty  as  the  resistance  tolerance  changes 
from  5%  to  0.1%  (@  R^^  tolerance  of  1%).  The  required  accuracy  in  absolute 
temperature  measurement  is  ±0.5°C.  Note  that  even  with  an  Rj^^  tolerance  of  1% 
and  an  R  tolerance  0.1%,  the  predicted  uncertainty  of  1.1“C  is  outside  system 


t 


requirements. 

The  25%  tolerance  on  reference-temperature  resistance  is  unacceptable. 

It  is  therefore  necessary  to  determine  experimentally  an  Rj^^  value  for  each 
thermistor.  This  implies  that  the  assumption  of  zero  uncertainty  in  the  value 
of  Tg  (Eq.  16)  is  no  longer  valid. 


Relative-Temperature  Uncertainty 


Because  of  the  highly  nonlinear  thermistor-resistance/temperature 
characteristics,  it  is  expected  that  there  will  be  uncertainty  in  AT 
measurements  as  well  as  absolute-temperature  measurements.  Beginning  with  the 
linear  approximation  (Eq.  3),  an  expression  can  be  obtained  for  AT 
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measurements  which  can  then  be  applied  to  Eq.  14  to  obtain  expressions  for 
relative-temperature  uncertainty . 

The  change  from  some  initial  temperature  to  some  later  temperature 
can  be  written  as 

r=r„-rj.  (3^) 


Solving  Bq.  3  for  T  and  combining  with  Eq.  31  we  obtain 

Ay- _ P  _ _ P  _ .  (32) 

(To In (i?i.y ■Rr.)+P)  (roln(J?r/i?r,)+P) 

Grouping  like  terms,  and  with  a  common  denominator,  we  produce  the  expression 


AT-p  I  +p)  (T,ln(i!,/Sr_)  .p) 

Applying  the  In  property;  ln(a)-ln(b)=ln(a/b)  to  Eq.  33  the  expression  is 


(33) 


Ar=- 


p  r|  ln(j;r/.Rr„) 


(34) 


( In  (  RtJRtJ  +  P  )  ( r«,  In  (  RtJ^t)  +  P  >  ' 

We  are  interested  in  determining  the  uncertainty  in  the  AT  measurements. 
In  Eq.  34  there  is  uncertainty  in  Ry,-,  R^„,  P,  and  T^. 

Applying  Eq.  34  to  Eq.  14,  the  uncertainty  in  AT  is  given  as 

11/2 


»ArH 


/aAr„\2jaAr„  ryaAr„ 


(35) 


where 


Tl  in-^ 

aAr_ 


ap 


r„in-^+p 

X\tp 


2P+T„ln 


Rj.^Rj.^ 


(36) 


T„ln.^+p)  (r„ln-yi+P 


dAT_ 

ar„ 


2 

Rrp  Rrp  Rji  Rtp 

3’o(2r„ln-^ln^+pln  ' 

To  To 

Rj.  \ 

r„ln-^+P 

■‘o  / 

Rj.  \ 

2 

(37) 
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Program  ERROR . RRS 

Equations  16-24  for  absolute-temperature  uncertainty  and  Eqs  35-40  for 
relative-temperature  uncertainty  were  incorporated  into  a  computer  program 
written  in  BASIC  entitled  ERROR. BAS.  The  program  gives  temperature 
uncertainty  in  the  initial  absolute  temperature  (presumably  averaged  over 
some  number  of  points  to  minimize  random  noise  effects),  as  well  as  the 
uncertainty  in  T^,  the  absolute-temperature  at  some  later  time,  and  the 
predicted  uncertainty  in  the  difference  of  the  two  (AT).  An  example  of  the 
results  of  an  analysis  is  given  in  Table  2.  Input  into  the  program  includes 
the  following: 
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P  =  3474  Kelvin, 

W.  =  104  Kelvin, 
t;  =  25'»C, 

W°  =  0.01®C  (assumed). 


r;”  =  10  kQ, 


WpT  varies  from  1%  to  25%, 

6^  °=  2.5  volts, 

VT  =0.01  volts, 

e®^T)  varies  with  temperature, 

W°  =  0.01%, 

R®°i'\0  kfi, 

W„ 


=  0  (assumed  that  a  close  tolerance  resistor 


is  available) . 

Additional  input  is  given  in  Table  2. 


Table  2.  Unnftrf.aintv  in  absolute  and  relative  temperature  with 
variation  in  reference  temperature  resistance  tolerance 
Temp  Input  Absolute  Temp _ Relative  Temperature - - 


h 

—n 

%[o 

% 

^in 

% 

44 

80 

25% 

7.44 

9.44 

2.35 

10% 

3.0 

4.28 

1.70 

5% 

1.64 

2.87 

1.58 

1% 

0.76 

2.23 

1.54 

34 

44 

25% 

6.8 

7.3 

0.68 

10% 

2.7 

3.0 

0.54 

5% 

1.4 

1.6 

0.52 

1% 

0.47 

0.76 

0.51 

A 

B 

C 

D 

E 

F 

Hlo 

2rIo 

%ri 

-RTn 

1.0 

0.0003 

1.8 

0.35 

1.03 

1.1 

0.0003 

0.71 

0.35 

1.03 

1.1 

0.0003 

0.35 

0.35 

1.03 

1.1 

0.0003 

0.07 

0.35 

1.03 

0.27 

0.00003 

0.45 

0.26 

0.35 

0.27 

0.00003 

0.18 

0.26 

0.35 

0.27 

0.00003 

0.09 

0.26 

0.35 

0.27 

0.00003 

0.02 

0.26 

0.35 

G 

H 

I 

J 

K 

The  U-terms  in  Table  2  (columns  G-K)  are  similar  to  those  defined  in  Eqs 
23-28.  They  are  the  individual  terras  in  the  relative-temperature-uncertainty 
equation  (Eq.  35)  included  here  in  order  to  determine  the  major  sources  of 
uncertainty  in  the  final  temperature  uncertainty  (columns  D-F) .  Two 
temperature  spans  are  considered  in  Table  2  in  order  to  observe  the  effects  on 
the  values  (column  F).  As  expected,  the  uncertainty  in  relative 
temperature  increases  with  increasing  temperature  interval.  For  each  case 
considered,  the  relative-temperature-uncertainty  cannot  be  considered 
negligible.  Note  also  that  the  error  in  P  (column  G)  becomes  comparatively 
substantial  as  the  tolerance  decreases  (for  example  in  rows  4  and  8). 
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Thermistor  Calibration  Techniques  for  Improved  Accuracy 


From  Tables  1  and  2  it  is  evident  that  the  25%  tolerance  is 
unacceptable  and  will  have  to  be  improved.  In  order  to  obtain  and 
values,  accurate  means  of  measuring  temperature  and  resistance  must  be 
determined.  A  number  of  different  schemes  are  considered  here,  including 
thermometer  temperature  measurements,  temperature  measurements  made  with  a 
precalibrated  thermistor,  resistance  measurements  using  a  digital  ohmmeter, 
and  resistance  values  obtained  through  voltage  measurements. 


Uncertainty  in  Rj^  Measurements 


In  order  to  obtain  R^^  at  reference-temperature  the  sensor  assembly 
must  be  brought  to  that  reference  temperature  and  the  resistance  recorded. 
Furthermore,  in  considering  the  sensor  thermal  time  constant  it  was  determined 
that  readings  should  be  taken  at  thermal  equilibrium.  It  is  the  practical 
problem  of  bringing  the  sensor  assembly  to,  say,  25®C  ±0.1®C  which  makes  this 
direct  approach  prohibitive. 

An  alternative  approach  is  to  take  a  reading  at  room  temperature  to 
obtain  T  and  R^,  then  designate  a  value  for  T^  and  use  Eq.  15  to  solve  for  Ry^,. 


R 


To 


(41) 


If  the  temperature  T  is  close  to  T^,  (for  example  T=22®C  and  Tq=  25®C)  the 
effects  of  nonlinear  P  on  Rj^  will  be  minimal.  However,  with  this  indirect 
approach  a  degree  of  uncertainty  is  inherent  in  the  resulting  Rj^,  due  to 
uncertainties  in  the  independent  variables  in  Eq.  41.  There  will  be 
uncertainty  in  T,  P,  and  R^.  In  order  to  quantify  the  overall  uncertainty  in 
the  calculated  R^^^  values,  an  uncertainty  analysis  is  performed.  Again,  the 
analysis  is  based  on  the  linear  approximation.  Applying  Eq.  41  to  Eq.  14  the 
uncertainty  equation  is  obtained: 
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where 


(42) 

(43) 

II 

1 

(44) 

(45) 

Equations  42-45  were  incorporated  into  a  computer  program  entitled 

RTNOT.BAS.  The  program  was  run  for  the  case  where  resistance  measurements  are 

made  using  a  digital  voltmeter  and  temperature  measurements  are  obtained  with 

a  National  Bureau  of  Standards  (NBS)  calibrated  thermometer  with  declinations 

of  0.1  C.  Input  into  the  program  includes: 

T  =  22®C  (ambient), 

W  =  0 . 1®C 

rJ  =  11,257  Q  (calculated  with  R^q=10,000  n), 

WpT  =  0.1%  of  reading  +  1  digit  =  20  n  or  0*2%'^, 

B  =  3474  K, 

Wg  =  104  K. 

RTNOT.BAS  gives  an  value  of  44  Q  or  0.4%  of  R^^ — a  substantial 

improvement  over  the  manufacturer's  tolerance  of  2500  fi  (25%  of  10,000  Q). 

The  individual  U-terms,  defined  as  the  product  of  sensitivity  Eqs.  43-45  and 
their  respective  tolerances,  give  Ug^=18,  Ug=35,  and  U^=20.  The  uncertainty  in 
P  is  now  the  leading  contributor  to  uncertainty  in  R^^. 

The  thermistor  sensitivity  in  the  temperature  range  of  interest  (35®C) 
is  about  235  Q/®C.  With  a  20  0  tolerance  on  the  ohmmeter  used  to  obtain  the 
resistance  measurements,  it  is  now  possible  to  read  resistance  variations 
relating  to  temperature  changes  of  0.1®C.  The  44  0  value  computed  for  the 
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case  above  relates  to  a  temperature  uncertainty  of  about  0.2®C.  As  a  check, 

ERROR. BAS  is  run  for  the  case  where 

T  =  T.  =  T  =  25®C 
R°  =  R  =  10,000  n 
e  =  2.5  volts 
B  =  3474  K 
Wp  =  104  K 
Wrto  =  44  £2 

Wg  =  10  £2  (0.1%  accuracy  assumed) 

W „  =  0.01  volts 
=  0.1®C. 

The  program  gives  a  0.25°C  uncertainty  in  Tj  and 

Thermistor  Calibration  Using  a  Precalibrated  Thermistor 

The  matching  of  thermal  time  constants  between  temperature  standard  and 
thermistor/sensor  is  an  important  consideration  in  thermistor  calibration.  To 
this  end  a  procedure  is  considered  which  utilizes  a  precalibrated  thermistor 
for  temperature  measurements. 

Initially,  one  or  a  number  of  thermistor/sensors  are  calibrated  as 
described  in  the  previous  section  using  an  ohmmeter  for  resistance 
measurements  and  a  MBS  calibrated  thermometer  for  temperature  measurements. 
Once  the  R^^  values  have  been  determined  for  the  reference  thermistor /sensors, 
they  can  be  used  to  measure  the  temperature  T  relating  to  the  resistance 
values  Rj  from  which  R^^  can  be  calculated  for  the  remaining  thermistor/ sensor 
assemblies.  Presumably  with  this  approach,  experimental  error  incurred  due 
to  mismatched  time  constants  will  be  reduced.  A  computer  program  was  written 
to  calibrate  thermistor/sensors  with  a  precalibrated  thermistor/sensor.  The 
program  RTNOTA.BAS  was  obtained  by  modifying  RTNOT.BAS  (see  Appendix).  It 
uses  the  reference-thermistor /sensor  R^  value  to  calculate  the  temperature, 
which  is  assumed  to  be  identical  for  both  the  reference—  and  test- 
thermistor/sensor  assemblies  (kept  in  close  proximity  during  testing).  This 
calculated  temperature  is  then  used  to  calculate  the  R^^j  value  of  the 
thermistor  being  calibrated. 
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Th©  pjTogjrain  usgs  Ecjs*  42*’45  to  calculat©  th.©  unc©]rtai.nty  in  the  Rjq 
measurements  so  that  for  each  run  the  expected  degree  of  accuracy  can  be 
monitored •  It  differs  from  program  RTNOT.BAS  in  that  it  uses  a  temperature 
uncertainty  (W^)  of  0.25®C  (as  determined  in  the  previous  section)  instead  of 
a  Wy  of  0.1®C  associated  with  direct  readings  from  the  NBS  calibrated 
thermometer.  A  major  concern  with  this  approach  is  that  the  0.25°C  Wj  value 
is  too  great  an  uncertainty  to  be  introduced  at  this  level  of  the  calibration 
procedure. 

An  analysis  was  performed  which  uses  ERROR. BAS  to  determine  the 
uncertainty  in  temperature  measurements  made  during  calibration  with  an  NBS 
calibrated  thermometer  as  compared  to  those  made  with  a  precalibrated 
thermistor/ sensor  assembly.  The  results  are  documented  in  Table  3. 
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Table  3.  Final  temperature  uncertainty  for  calibration 
temperature  measurements  obtained  with  a  thermometer  compared 
to  those  made  with  a  thermistor. 


Thermometer 

Thermistor 

A 

B 

C 

D 

E 

F 

Input 

T 

22 

25 

25 

22 

25 

25 

T? 

25 

25 

25 

25 

25 

25 

t’ 

25 

44 

80 

25 

44 

80 

r" 

10,000 

10,000 

10,000 

10,000 

10,000 

10,000 

RT^ 

11,258 

10,000 

10,000 

11,258 

10,000 

10,000 

o 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

w! 

104 

104 

104 

104 

104 

104 

WpTn 

20 

40 

40 

20 

40 

40 

Wp 

10 

10 

10 

10 

10 

10 

wL 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

wr 

0.05 

0.24 

0.24 

0.24 

0.34 

0.34 

Output 


0.24 

0.34 

0.33 

0.34 

0.41 

0.41 

Wn 

0.24 

0.76 

2.23 

0.34 

0.81 

2.26 

”at 

0.29 

0.64 

1.91 

0.30 

0.64 

1.91 

U-Terms 

Ud 

0 

0.5 

1.6 

0 

0.5 

1.6 

0 

0.001 

0.009 

0 

0.001 

0.01 

OpT 

0 

0.01 

0.04 

0 

0.01 

0.04 

0.21 

0.21 

0.21 

0.21 

0.21 

0.21 

'^RTn 

0.21 

0.35 

1.1 

0.21 

0.35 

1.1 

Note:  Beta  held  constant  throughout  at  3474  K,  and  input  values  in 
columns  A  and  D  were  calculated  using  the  linear  approximation. 


Columns  A,  B,  and  C  are  for  the  case  of  thermometer  temperature 
measurements  and  columns  E,  and  F  are  for  thermistor  temperature 
measurements.  Column  A  depicts  the  case  where  the  value  (@  25®C)  of  a 
thermistor  is  being  determined  from  resistance/temperature  data  obtained  at 
room  temperature  (22®C).  The  input  temperature  uncertainty  (W^)  in  column  A 
is  O.OS^C.  (Here  it  is  assumed  that  the  0.1®C  declinations  on  the  thermometer 
can  be  read  within  ±0.05°C.)  The  value  of  20  0  assumes  resistance 

measurements  made  with  a  digital  ohmmeter  and  the  Wp  value  of  10  Q  assumes  a 
0.1%  resistor  tolerance.  The  output  shows  an  absolute-temperature  uncertainty 
of  0.24®C.  This  value  is  used  as  input  for  column  B.  It  means  that  the 
uncertainty  in  the  T^  value  of  25°C  is  0.24°C.  Note  also  that  the  reference 
temperature  resistance  tolerance  {W|^i-q)  has  increased  from  20  Q  in  column  A  to 
40  Q  in  column  B;  this,  too,  is  due  to  solving  for  the  thermistor  R^^  value 
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from  data  taken  at  22®C  rather  than  obtaining  the  value  directly  at  25'>C. 

The  output  for  column  B  shows  an  increasing  temperature  uncertainty  from 
that  of  column  A.  Note  that  the  value  relates  to  the  input  value  of  44°C. 
The  only  difference  between  columns  B  and  C  is  the  increased  temperature  span 
to  80®C.  In  all  three  columns  the  relative-temperature  uncertainty  is 
substantial.  The  U-Terms  show  P  to  be  the  major  contributor  to  final 
temperature  measurement  uncertainty. 

In  examining  the  values  in  columns  D,  E,  and  F,  which  relate  to  the  case 
of  thermistor  temperature  calibration  with  a  precalibrated  thermistor,  it  can 
be  seen  that  the  overall  temperature  uncertainty  shown  in  the  output  is 
greater  than  in  columns  A,  B,  and  C.  In  column  D  the  temperature  uncertainty 
input  (W^)  is  given  as  0.24‘'C.  This,  as  in  column  A,  relates  to  the  case 
where  the  (@  25°C)  value  of  the  thermistor  is  being  determined  indirectly 
from  temperature/resistance  data  obtained  at  ambient  temperature  (22‘*C). 
Because  the  thermistor  depicted  in  column  D  is  being  calibrated  with  a 
precalibrated  thermistor/sensor  assembly,  it  is  exactly  the  same  case  as  that 
depicted  in  column  B  with  the  exception  that  the  value  is  20  n 

(measurement  with  ohmmeter).  The  thermistor  has  not  yet  acquired  the 
uncertainty  which  is  introduced  by  indirectly  obtaining  values  from  data 
obtained  at  22 ®C.  It  is  not  until  column  E  that  this  uncertainty  is 
introduced . 

Again,  the  only  difference  between  the  input  values  in  columns  E  and  F 
is  the  T^  values  of  44®C  and  80®C,  respectively.  The  output  temperature 
errors  shown  in  columns  B,  C,  E,  and  F  are  simply  too  great  to  introduce 
during  thermistor  calibration.  The  use  of  a  precalibrated  thermistor  for 
thermistor  calibration  is  not  practical  with  the  temperature  and  resistance 
measuring  techniques  considered  here.  The  U-term  values  given  in  columns  E 
and  F  show  P  as  the  major  source  of  final  temperature  uncertainty. 
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Experinental  Determination  of  Beta  for  Improved  Accuracy 


It  can  be  seen  in  Tables  3  and  4  that  the  greatest  error  incurred  in 
temperature  measurements  is  due  to  uncertainty  in  P*  Thermometries  reports 
that  p  varies  with  a  production  lot  by  3%  and  from  lot  to  lot  by  up  to  5%.  We 
have  been  assuming  a  3%  variation.  At  this  point  this  represents  the  major 
limit  in  our  ability  to  accurately  calibrate  thermistor/sensor  assemblies.  It 
has  been  reported  that  P  should  be  considered  unique  to  each  thermistor  and 
should  be  experimentally  determined^^.  Whereas  it  is  necessary  to 
experimentally  determine  through  resistance/temperature  measurements  at 
one  or  a  number  of  points,  it  is  possible,  at  the  same  time,  to  use  this  data 
experimentally  to  determine  p. 

In  order  to  establish  the  accuracy  with  which  P  can  be  found 
experimentally,  an  uncertainty  analysis  is  performed  based  on  the  linear 
approximation.  Solving  Eg.  3  for  P  we  obtain 


T  T 

B=  q. 


(46) 


Uncertainty  in  experimentally  determined  P  is  expected  as  a  result  of 
uncertainty  in  the  independent  variables  R^,  T  and  T^.  The  uncertainty 

equation  (Eq.  14)  becomes 


where 


(47) 


(48) 

dR^  Rj.(T^-T) 


ap  _  -TT„ 

dR^;  ’ 


(49) 
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(50) 


dT 


dT^ 


(r„-D 


(51) 


These  expressions  were  incorporated  into  a  computer  program  ER_BETA.BAS, 
and  a  study  was  performed  using  this  program  to  determine  experimental  P 
uncertainty.  The  results  are  given  in  Table  4. 

From  Eq.  46  it  is  evident  that  in  order  to  calculate  P  using  the  linear 
approximation,  two  temperature/resistance  data  points  must  be  experimentally 
obtained.  If  P  is  temperature  dependent,  then  it  is  important  to  determine 
the  temperature  range  over  which  we  are  concerned.  It  has  been  concluded  that 
this  temperature  range  is  between  22 “C  and  50°C.  Determining  p,  the  slope  of 
the  In(Rj)  versus  1/T  plot,  between  22®C  and  50®C  will  give  us  our  most 
accurate  temperature  measurements  at  a  point  halfway  between  those  two  points, 
or  in  the  36®C  vicinity.  Similarly,  determining  P  between  36  and  50®C  we 
determine  the  slope  in  the  43 °C  region  and  obtain  our  most  accurate  readings 
in  that  region.  With  this  in  mind,  p  uncertainty  was  determined  for  a  number 
of  cases.  In  columns  A  and  B  of  Table  4  the  reference-temperature  is  taken  to 
be  36»C  with  an  additional  point  at  22°C.  Note  also  the  difference  in 
temperature  measurement  uncertainty  in  the  input  variables  W^^  and  W^.  In 
column  A  it  is  assumed  that  the  NBS  calibrated  thermometer  can  be  read  to 
±0.1®C.  In  column  B  a  ±0.05°C  readability  is  assumed.  In  the  former  case  the 
U-Terms  show  temperature  measurement  uncertainty  is  the  primary  source  of  P 
uncertainty.  In  the  latter  case  it  is  our  ability  to  measure  resistance 
accurately  (again,  it  is  assumed  that  resistance  measurements  are  obtained 
with  a  digital  ohmmeter). 
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A  B  C  D  E  F 

36  36  22  36  50  36 

22  22  36  50  36  22 

0.1  0.05  0.05  0.05  0.05  0.05 

0.1  0.05  0.05  0.05  0.05  0.05 

20  20  20  20  20  2 

20  20  20  20  20  2 


Output 

We 

42 

29 

29 

45 

45 

17.7 

%^g 

1.2 

0.8 

0.8 

1.3 

1.3 

0.51 

U-Terms 

12 

12 

20 

35 

22 

1.2 

Upl. 

20 

20 

12 

22 

35 

2.0 

26 

13 

12 

12 

13 

13 

24 

12 

13 

13 

12 

12 

Note:  Beta  held  constant  throughout  at  3474  K. 


Alternating  reference  temperatures  has  no  effect  on  P  uncertainty  as  can 
be  seen  in  the  output  Wg  values  in  columns  B  and  C.  Similarly,  in  columns  D 
and  B  temperatures  of  36‘’C  and  50»e  are  alternated  as  reference  temperature 
with  the  same  result.  Beta  uncertainty  does  increase  at  elevated  temperature 
intervals  as  seen  by  the  increase  in  Wg.  Output  %Wg  terms  should  be  compared 
to  the  manufacturer's  p  tolerance  of  3%.  In  column  F,  resistance  measurement 
uncertainty  is  decreased  in  order  to  determine  the  effect  on  our  ability  to 
measure  p.  The  improvement  is  substantial  with  a  60%  decrease  in  P 
uncertainty,  and  a  95%  decrease  in  U-Terms;  Ugj  and  Overall  system 

accuracy  is  greatly  improved  with  more  accurate  resistance  measurements. 

The  question  now  arises  as  to  how  this  increased  beta  precision  affects 
the  overall  temperature  measurement  uncertainty.  A  study  was  performed  using 
program  ERROR. BAS  for  this  purpose.  The  program  was  modified  to  include  in 
the  output  the  U-Terms  from  the  absolute-uncertainty  expressions  as  well  as 
those  from  the  relative-temperature  uncertainty  expressions.  The  results  of 
this  study  are  documented  in  Table  5. 
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Table  5.  Effects  of  increased  B  accuracy  on 


Input 

!i 

k 


RTO 

To 


Output 


W. 


W 

W' 


Tn 


U-Terms  for 


Up 

u 


u 


RTo 

RTi 


RTn 


U^Terms  for  W. 


B 

To 

RTo 

RTi 


To 

RTo 

RTn 


A 

B 

C 

D 

36 

36 

36 

36 

22 

30 

22 

30 

50 

45 

50 

45 

6603 

6603 

6603 

6603 

2.5 

2.5 

2.5 

2.5 

35 

35 

17 

17 

20 

20 

2 

2 

0.05 

0.05 

0.05 

0.05 

0.25 

0.26 

0.20 

0.24 

0.45 

0.39 

0.43 

0.37 

0.50 

0.44 

0.47 

0.43 

CM 

-T 

0.1 

0.1 

0.05 

0.00003 

0.00004 

0.00003 

0.00004 

0.015 

0.008 

0.0015 

0.0008 

0.19 

0.23 

0.19 

0.23 

0.42 

0.36 

0.42 

0.35 

-S?Ti 

I  ^  A 

-0.13 

0.06 

0.07 

0.03 

0.05 

0.05 

0.05 

0.05 

0.08 

0.08 

0.008 

0.008 

0.19 

0.23 

0.19 

0.23 

^.15 

0.09 

0.07 

0.05 

0.05 

0.05 

0.05 

0.05 

0.09 

0.09 

0.009 

0.009 

0.41 

0.36 

0.42 

0.36 

Note;  Input  Rj^  calculated  using  linear  approximation, 
constant  throughout  at  3474  K. 


Beta  held 


In  columns  A  and  B  it  is  assumed  that  resistance  is  measured  with  an 
ohmmeter.  Columns  C  and  D  assume  a  more  accurate  resistance  measurement 
capability,  which  manifests  itself  in  lower  input  Wg  and  values.  The 

temperature  interval  for  columns  A  and  C  is  22 °C  to  50°C.  For  columns  B  and  D 
it  is  30®C  to  45®C.  In  each  case  the  reference  temperature  is  36®C. 
Temperature  measurements  are  assumed  to  be  made  with  a  thermometer.  The 
higher  P  accuracy  is  evidenced  by  a  decrease  in  absolute  and  relative 
temperature  uncertainty  now  shown  to  be  0.50®C  and  below.  U~terms  show  the 
primary  source  of  temperature  measurement  uncertainty  to  be  in  the  resistance 
measurements  rather  than  the  P  uncertainty. 
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If,  in  calibrating  the  thermistor/ sensor  assemblies, 
temperatqre/resistance  measurements  are  obtained  at  three  temperatures,  Eq,  46 
can  be  solved  for  3  different  cases: 

@  36°C 

T,  e  22°C 

2-  Rto  0  36°C 

T,  Ry  @  50<*C 

3.  T^,  Rto  0  22»C 

T,  Ry  @  50°C. 

Where  the  subscripts  L,  M,  and  H  denote  low,  medium,  and  high  temperatures, 
respectively.  Additionally,  a  Pgverage  determined  using  the 

expression 


^average 


3 


(52) 


This  value  is  expected  to  be  very  close  to  the  value  of  8^,). 

A  computer  program  was  developed  to  calculate  these  four  P  values.  The 
program,  EX_BETA.BAS,  also  uses  Eqs.  47-51  to  determine  the  uncertainty  in 
each  P  value.  This  ability  has  been  included  in  the  program  so  that 
variations  in  the  four  P  values  can  be  compared  to  the  predicted  uncertainty 
in  each  value.  If  P  value  variations  are  consistently  outside  the  predicted 
uncertainties,  this  can  be  taken  as  evidence  of  a  temperature-dependent  P 
which,  if  a  high  degree  of  accuracy  is  needed,  can  be  accounted  for  by  some 
polynomial  P  expression. 


Determination  of  Thermistor  Resistance  Using  Voltage  Measurements 


It  has  been  shown  (Table  5)  that  the  primary  source  of  temperature 
uncertainty  is  in  the  resistance  measurements  when  the  temperature/resistance 
values  are  obtained  with  a  NBS-calibrated  thermometer  and  digital  ohmmeter. 
The  ±20  Q  meter  accuracy  causes  uncertainty  in  both  the  R^^  and  P  values. 
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Another  approach  is  considered  which  utilizes  voltage 
measurements  obtained  using  Labtech  Notebook.  Resistance  values 
are  then  computed  using  Eq.  11: 


(11) 


In  addition,  the  source  voltage  can  be  recorded  simultaneously.  Source- 
voltage  drift  is  a  problem  that,  up  to  now,  has  not  been  addressed.  Using 
Eqs.  21—24  the  uncertainty  in  the  Rj.  measurements  for  this  case  was  predicted. 


Using  the  values 


e  =2.57  volts  (obtained  experimentally), 

W*  =  2.57x10"^  volts  (0.01%  tolerance), 
e®^T)  =  1.32  volts  (obtained  experimentally), 

W®  =  1.32x10'^  volts  (0.01%  tolerance), 

R®°i^^,998.8  n, 

Wp  =  0.2  Q  (1%  tolerance  resistor  calibrated  at  Natick 
"  RD  &  E  electronics  lab  to  a  tolerance  of  0.002%). 


Eq.  14  gives  a  value  of  2.8  Q  (0.03  %  tolerance);  a  substantial 
improvement  over  the  20  a  tolerance  obtained  with  ohmmeter  resistance 
measurements.  When  the  new  Rj  tolerance  is  taken  into  account,  the  final 
temperature  measurement  accuracy  is  greatly  increased.  Prom  ERROR. BAS  we 
obtain  the  following  numbers; 


Table  6.  Temperature  uncertainties  for  case  of  resistance  values 
obtained  through  voltage  measurements  >_ 


Input 

A 

B 

p 

3474 

3474 

T„ 

36 

36 

? 

Tn 

25 

25 

50 

80 

Rt 

To 

6603 

6603 

2.57 

2.57 

Wb 

”to 

18 

18 

2.8 

2.8 

0.05 

0.05 

Output 

^AT 

0.07 

0.07 

0.09 

0.27 

0.096 

0.23 

Two  temperature  ranges  are  considered  in  Table  6;  in  column  A  the 
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temperature  span  is  25®C  -  50®C  with  the  maximum  output  uncertainty  (0.096°C) 
appearing  in  the  relative-temperature  term,  while  in  column  B  the  temperature 
span  is  25°C  -  80°C  with  the  maximum  output  uncertainty  (0.27®C)  appearing  in 
the  high-temperature  term.  The  input  value  of  0-05°C  assumes  temperature 
measurements  made  with  a  thermometer.  The  U-terms  (not  included)  show  the 
primary  source  of  error  to  be  in  |5  due  to  temperature  measurements  of  ±0.05®C 
accuracy. 

It  is  important  to  remember  that  the  uncertainty  analysis  is  based  upon 
the  constant  P  linear  approximation.  For  the  25®C  -  80°C  temperature  range  in 
column  B,  the  actual  error  in  an  80°C  measurement  could  be  much  greater 
(±0.4®C)  than  that  shown  in  due  to  variations  in  p.  For  our  temperature 
range  (20®C  -  50°C),  however,  P  variations  are  expected  to  introduce  an 
additional  error  of  ±  0.1°C.  When  this  is  added  to  the  output  in  column  A  we 
are  Still  inside  accuracy  requirements,  and  the  determination  of  thermistor 
resistance  through  voltage  measurements  represents  an  acceptable  thermistor 
calibration  approach. 

RTNOX.BAS  Output  for  Voltage  Rj.  Measurements 

In  the  section  entitled  Uncertainty  in  Measurements,  a  procedure  was 
investigated  whereby  the  value  of  a  thermistor  is  calculated  using  the 
linear  approximation  from  temperature/resistance  values  recorded  at  a 
temperature  close  to,  but  not  equal  to,  T^.  It  was  found  that  this  procedure 
introduces  additional  uncertainty  into  the  value.  RTNOT.BAS  gave  a 
value  of  44  ohms  for  this  case. 

With  more  accurate  resistance  values  obtained  through  voltage 
measurements,  this  procedure  might  now  be  feasible.  An  experimental  run  is 
performed  where  temperature/ resistance  data  are  taken  at  22.20®C.  Using  these 
data  Rt  is  calculated  for  T„  =  25°C.  RTNOT.BAS  is  run  again  with  the 
following  results: 
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Input 

T  (ambient)  =  22.20®C 
W,  =  0.05“C  (thermometer) 

P  =  3474  K  (approximation) 

R®  (ambient)  =  11,415.5  ohms  (experimental) 
Wgj  (ambient)  =  3  ohms. 


Output 

R,  =  10,222  Q  0  25®C 

Tolerance  =  2.22%  (relative  to  R^^  =  10,000  n,  as  purchased) 
WoT„  =  21.3  n 

=  0.21%  (relative  to  R^^,  =  10,222  0) 


U-Terms 
Urt  =  2.7 
Or  =  5.9 

Uj  =  20 


The  output  now  shows  a  Wr^,,  value  of  21.3  Q,  a  100%  improvement  over  the 
previously  determined  44  Q  case,  but  far  from  the  3  Q  uncertainty  obtained 
when  the  Rj^  value  is  acquired  at  the  exact  temperature  to  which  it  relates. 


Temperature  Correction  Terms 

Temperature  measurement  error  is  introduced  in  the  form  of  a  bias  when 
power  dissipated  in  a  thermistor  causes  it  to  self-heat.  This  section 
describes  the  formulation  of  a  corrective  term,  which  eliminates  this  bias. 

In  addition,  an  expression  is  presented  for  an  emergent  stem  correction  which 
is  required  when  a  thermometer  is  used  as  the  calibration  standard. 


Experimental  Investigation  of  Thermistor  Self-Heating 


Power  is  dissipated  in  a  thermistor  in  the  form  of  heat  when  connected 
to  an  electrical  circuit.  This  causes  the  thermistor  temperature  to  rise 
above  the  ambient,  introducing  an  error  between  the  actual  ambient  temperature 
and  the  temperature  of  the  thermistor .  An  energy  balance  can  be  written  as 

^^lost  j.  ^Oabsorbed  (53) 

dt  dt  dt 

This  states  that  the  rate  at  which  energy  is  supplied  to  a  thermistor  is  equal 
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to  the  rate  at  which  it  is  lost  to  its  surroundings  plus  the  rate  at  which  it 
is  absorbed.  The  rate  at  which  energy  is  supplied  is  equal  to  the  power 
dissipated  in  the  thermistor 


te„(r)]  J  =  Rj.= 


[e„(D]^ 


(54) 


where  e  (T)  is  the  voltage  across,  and  Rj,  the  resistance  in  the  thermistor. 
The  rate  at  which  heat  is  lost  is  proportional  to  the  thermistor  temperature 
rise; 

(r-Ta)  . 

dt  ^ 

Where  6  is  the  dissipation  constant,  T  is  the  thermistor  temperature,  and  T^ 
is  the  ambient  temperature.  Power  is  absorbed  by  a  thermistor  according  to 
the  expression 

(56) 

~dr  dt’ 


where  Cp  is  the  specific  heat  of  the  thermistor  material  and  m  is  its  mass. 
Combining  Eqs.  53-56  we  obtain  the  expression 


dt  Rt  '  ^  p  dt 


(57) 


After  a  relatively  short  period  of  time  (t  >  thermal  time  constant  = 
0.21  s),  the  thermistor  is  at  thermal  equilibrium  and  the  absorption  term 
drops  out  of  Eq.  57  leaving  the  expression; 

P=6(T-Tg).  (58) 


In  order  to  quantify  the  thermistor  temperature  rise  due  to  self¬ 
heating,  the  dissipation  constant  6  must  be  determined.  If  two  power  readings 
are  taken  we  obtain  two  expressions 

Pi=6(Ti-Ta)  (58) 


anci 
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P2=1>{T^-TJ  . 


(60) 


In  order  to  eliminate  T^,  Eq.  60  is  subtracted  from  Eq.  59,  which  is  valid 
providing  the  ambient  temperature  does  not  change.  The  expression  becomes 

or 

Power  terms  can  be  calculated  using  the  relation 

Rj.  ' 

where,  for  a  linear  voltage  divider 

where  e^  is  the  source  voltage  and  R  is  the  resistance  of  the  voltage  divider 
resistor.  If  the  6  value  is  now  inserted  into  Eq.  58,  the  temperature  rise 
corresponding  to  a  given  power  can  be  computed; 

^=Tsir=(T-Tj  .  (65) 

Table  7  gives  the  results  of  an  experimental  test  run  to  determine  6  for 
thermistor/ sensor  assembly  Tjg. 


(61) 


(62) 


(63) 


(64) 


Temp 

21.85 


Table  7.  Experimental  determination  of  dissipation 
constant  6  at  21.85°C. 


1.094 

0.^89 

11^00 

T 

21.59 

P 

2.96E-5 

PZ6 

0.087 

T  -  P/5 
21.50 

2.257 

1.211 

11588 

21.86 

1.27E-4 

0.375 

21.49 

3.781 

2.005 

11308 

22,54 

3.56E-4 

1.052 

21.50 

4.327 

2.281 

11175 

22.88 

4.66E-4 

1.380 

21.50 

The  temperature  shown  in  the  first  column  is  that  of  the  thermometer  and 
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was  held  constant  fpr  each  of  the  four  runs-  In  order  to  ensure  a  constant 
Tg,  the  system  was  brought  to  thermal  equilibrium  and  source-voltage  was 
turned  on  only  while  data  were  being  acquired-  The  source-voltage  was 
varied  and  the  voltage  across  the  thermistor  was  calculated  for  each  case 
using  Eq.  64.  The  variation  in  Ry  values  in  Table  7  shows  the  error 
introduced  by  power  dissipation-  This  is  also  evidenced  in  the  rising 
temperatures jr  which  were  calculated  using  the  linear  approximation  and 
previously  determined  P  and  values-  Power  was  calculated  using  Eq-  63. 
The  P/5  entries  are  from  Eq-  65  and  have  units  of  degrees  Celsius-  These  are 
the  corrective  values-  When  the  corrective  term  is  subtracted  from  the 
calculated  temperature,  the  final  temperature  is  approximately  the  same  for 
each  run.  These  corrected  temperatures  should  be  the  same  as  the  thermometer 
temperature  (in  this  case  they  are  not  because  the  correction  was  not 
implemented  in  the  original  determination  of  P  and  R^^) . 

The  value  of  the  dissipation  constant  used  above  was  determined  by 
applying  Eq-  62  first  to  the  data  obtained  in  runs  1  and  2,  then  to  the  data 
obtained  in  runs  1  and  3,  etc.  until  all  possible  combinations  were  exhausted. 
The  value  of  6  determined  for  each  combination  is  shown  in  Table  8- 


Table  8-  Dissipation  constant  values  determined 


Combination 

Temp  =  21.85«>C  1-2 

1-3 

1- 4 

2- 3 

2- 4 

3- 4 


using  data  given  in  Table  7- 

Dissioation  Constant  (W/°C) 
3.58x10'^ 

3.39x10*^ 

3.37x107 
3.33x10’^ 

3.32x10'^ 

3.30x10'^ 

3, 38x10'^  W/<»C 
1.02x10'5  W/®C 


Average  value  used  in  Table  7  = 
Standard  deviation  = 


The  study  shown  in  Tables  7  and  8  was  repeated  at  a  higher  temperature 
to  determine  the  degree  to  which  the  dissipation  constant  is  temperature 
dependent-  The  results  are  given  in  Tables  9  and  10. 
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Table  9.  Experimental  determination  of  dissipation 
constant  5  at  49.15”C. 


Temp  Run 

49.15  1 

2 
3 


0.I21  o.ieo 

2.676  0.839 

4.549  1.400 


R  T 

46^7.7  48.31 

4563.3  48.74 

4452.7  49.49 


£  ZLA 

1.46E-5  0.041 
1.54B-4  0.434 
4.40E-4  1.239 


T  -  P/S 
48.27 
48.31 
48.25 


Table  10.  Dissipation  constant  values  determined  using  data 


Combination 

Temp  =  49.15°C  1-2 

1- 3 

2- 3 


given  in  Table  9> 

Dissipation  Constant  W/°C 
3.24x10’^ 

3.60x107 

3.80x10'^ 


Average  value  used  in  Table  9  =  3.55x10  ^  W/*C 
Standard  deviation  =  2.84x10'^  W/®C 


The  relatively  small  (5%)  change  in  5  from  Table  8  to  Table  10  is  an 
indication  that  the  value  of  5  is  not  greatly  affected  by  changes  in 
temperature  (though  it  would  rec[uire  a  larger  number  of  samples  to  determine  a 
statistically  sound  basis  for  that  assumption).  The  amount  of  power 
dissipated  in  a  thermistor/ sensor,  however,  is  strongly  temperature  dependent. 


Fig.  2.  Self-heating  temperature  correction  curve 
for  thermistor/ sensor  Tjg. 
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Figure  2  is  a  plot  of  corrective  temperature  (P/5)  versus  absolute- 
temperature  for  thermistor/ sensor  Tjg.  The  plot  shows  the  error  incurred  if 
self-heating  is  not  taken  into  account.  Upon  subtraction  of  this  curve,  the 
data  will  relate  to  the  desired  ambient  temperature. 

Figure  2  was  generated  using  the  linear  approximation  (Eg.  3).  Power 
dissipated  was  determined  using  Eqs.  63  and  64.  The  correct ive-temperature 
term  was  calculated  using  Eg.  65.  In  addition,  the  following  values  were 
used: 

=  3368  K  (experimental), 

=  7003  ohms  (experimental) 

@  Tq  =  36°C, 

5  =  3.3x10'^  Watts/^C, 

Gg  =  2.5  V. 

The  dissipation  constant  has  been  previously  given  as  0.045  mW/®C  in 
still  air  and  0.23  mW/°C  in  the  water  plunge.  The  nominal  value  computed  here 
of  0.33  mW/^C  represents  a  decrease  in  temperature  rise  due  to  self-heating 
when  the  thermistors  are  mounted  in  sensors. 

In  Tables  8  and  10  standard  deviations  of  0.01  mW/®C  and  0.028  mW/°C, 
respectively,  are  given.  In  order  to  determine  the  significance  of  these 
variations,  an  uncertainty  analysis  is  performed.  Combining  Eq.  14  with  Eg. 
62,  an  expression  is  obtained  for  the  uncertainty  in  6: 


where 


(66) 


(67) 
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(70) 


The  uncertainty  in  T^  and  T2  is  known  to  be  approximately  0.1®C.  The 
values  of  Wp,  and  Wpj  can  be  obtained  by  combining  Eq.  14  with  Eq.  63: 


dP 

3e„(T) 


(71) 


where 


dp  _2[e„(D] 

ae„(D  ■  J?j, 


and 


[e„(r)]^ 
dRr  (Rt)^ 


The  uncertainty  in  thermistor  voltage  (VI is  unknown  but  can  be 
computed  by  combining  Eq.  14  with  Eq.  64.  The  expression  is 


where 
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Using  the  experimental  data  of  runs  1  and  4  in  Table  7  and  Eqs.  67-77, 
the  uncertainty  in  6  is  computed  by  Eq.  66  to  be  0.023  mW/®C.  This  is  quite 
close  to  the  computed  experimental  standard  deviations  of  0.01  mW/®C  and 
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0.028  mW/®C,  and  is  taken  as  an  indication  that  this  is  a  reliable  technique 
by  which  to  determine  the  dissipation  constant  5.  In  terms  of  temperature, 
this  0.023  mW/®C  uncertainty  in  the  dissipation  constant  corresponds  to 
±0.04®C  at  25 and  a  source-voltage  of  2.5  V. 

Because  the  dissipation  constant  is  highly  dependent  on  heat  transfer 
mechanisms,  which  can  vary  from  one  thermistor  to  another,  it  is  necessary  to 
determine  a  unique  dissipation  constant  for  each  thermistor/ sensor .  Program 
EX  BETA. BAS  was  modified  to  include  the  calculation  of  6.  The  new  program, 

DC  BETA. BAS,  requires  that  an  additional  run  be  made  at  low  temperature,  which 
will  provide  the  information  necessary  to  solve  Eqs.  62-65.  This  run  is 
p©r formed  at  a  lower  source— voltage  (power)  than  that  used  to  determine  P  and 

^0- 

As  can  be  seen  from  Eq.  62,  in  order  to  calculate  5,  the  temperatures  T-| 
and  T2  relating  to  powers  P-j  and  P2  must  be  calculated  from 
temperature/ resistance  data  obtained  in  the  same  run. 


(62) 


In  order  to  determine  these  temperatures,  P  must  be  known  for  that  particular 
thermistor/ sensor.  Beta  cannot  be  accurately  determined  without  the 
self-heating  correction. 

For  this  reason  an  iterative  technique  is  implemented  whereby  a  value  of 
6  is  assumed,  self-heating  is  computed  based  on  that  assumption,  temperature 
is  corrected,  P  is  determined,  corrected  temperatures  are  calculated,  which 
are  then  used  to  determine  a  new  5,  etc.  The  process  continues  until  the 
difference  between  any  two  consecutive  values  of  5  is  less  than  1x10"^  W/°C. 
This  value  is  approximately  an  order  of  magnitude  less  than  the  predicted  6 
uncertainty. 

With  the  value  of  5  known  for  each  thermistor/sensor  the  self-heating 
correction  can  be  implemented,  depending  on  the  required  accuracy,  either  by 
assuming  an  average  5  value  for  all  thermistor /sensor  assemblies,  or  by 
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allowing  each  to  retain  its  own  unique  value.  In  either  case,  a  single 
expression,  obtained  by  combining  Eqs  63-65,  can  be  used  to  obtain  the 
'actual'  ambient  temperature; 


(78) 


The  alternative  to  implementing  corrections  for  thermistor  self-heating 


is  to  operate  at  a  source  voltage  for  which  thermistor  self-heating  can  be 
considered  negligible.  From  Fig.  2  it  can  be  seen  that  the  maximum  error  due 
to  thermistor  self-heating  occurs  at  about  25‘'C.  Fig.  3  shows  a  curve  which 
was  generated  using  the  second  term  on  the  right  of  Eq.  78; 


Rji 

■y  ■ 


(79) 


An  Rj  value  relating  to  a  temperature  of  25®C  (10,000  Q)  was  held 
constant.  A  value  of  3.3x10’^  W/®C  was  used  for  the  dissipation  constant,  an 


Fig  3.  Temperature  error  versus  source  voltage  for 
thermistor  self-heating  at  25®C  and  6  -  3.3x10"^  W/®C. 

R  value  of  10,000  Q  was  used,  and  the  source  voltage  was  allowed  to  vary  from 
0  volts  to  2.5  volts.  The  curve  shows  that  a  source  voltage  of  1.2  V  or  less 
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would  keep  the  maximum  error  inside  that  of  experimental  uncertainty.  For  a 
source-voltage  in  the  vicinity  of  1  V,  experimental  data  have  shown  thermistor 
voltages  ranging  from  about  0.2  V  to  1  V  for  the  temperature  range  of 
interest . 

Thermometer  C!Qiergent*^Stem  Correction 

The  NBS  calibrated  thermometer  used  as  the  calibration  standard  is  a 
total-immersion  thermometer.  It  is  not,  however,  being  used  as  such.  A 
correction  must  be  made  to  compensate  for  that  portion  of  the  thermometer 

exposed  to  room  temperature  when  readings  are  obtained  above  room  temperature. 

.  •  •  14 

The  commonly  used  equation  for  this  correction  is 

T^^=Kn{T-t),  (80) 

where  K  is  the  differential  expansion  coefficient  of  mercury,  n  is  the  number 
of  degrees  emergent  from  the  bath,  T  is  the  temperature  of  the  bath,  and  t  is 
the  mean  temperature  of  the  emergent -stem.  Here,  the  values  are 
K  =  0.00016, 

n  =  lll^C, 

T  =  36®C  and  BO^C, 

t  :  determined  experimentally. 

The  value  of  t  is  determined  by  fastening  another  thermometer  to  the  stem  of 
the  standard  at  approximately  the  halfway  point.  With  t  =  25 °C  the  emergent- 
stem  correction  is  0.20®C  @  and  0.44“C  @  50®C.  Both  of  these  values  are 

outside  our  predicted  temperature  measurement  uncertainty  and,  therefore,  must 
be  taken  into  consideration  during  thermistor/sensor  calibration. 

Suumaxy  of  Uncertainty  Analysis 

It  has  been  shown  that  in  order  to  achieve  the  required  temperature 
measurement  precision  of  ±0.5®C,  a  thermistor  calibration  procedure  must  be 
implemented.  Furthermore,  it  is  possible  to  attain  the  required  tolerances 
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with  equipment  that  is  currently  available  at  Natick.  The  procedure  requires 
the  acquisition  of  four  temperature/resistance  data  points  at  three  different 
temperatures.  The  temperatures  correspond  to  the  midpoint  and  extremes  of  the 
temperature  range  of  interest,  which  was  designated  to  be  22 "C  -  BO^C 
(midpoint  at  36“C) . 

The  dissipation  constant  is  determined  at  22 ®C  by  taking  readings  at  two 
power  settings  relating  to  source-voltage  inputs  of  0.5  V  and  2.5  V. 

Obtaining  the  dissipation  constant  at  this  temperature  ensures  the  greatest 
degree  of  accuracy  because  of  the  high  thermistor  sensitivity  to  self-heating 
at  this  temperature  (greatest  change  in  dissipated  power  per  unit  change  in 
source-voltage) .  The  uncertainty  analysis  predicts  variations  in  6  to  be  in 
the  vicinity  of  a  standard  deviation  of  0.025  mW/®C. 

Reference-temperature  resistance  (Rj,,)  relating  to  a  reference 
temperature  (T^^)  of  36°C  is  calculated  from  temperature/voltage  data  taken 
close  to,  but,  in  general,  not  equal  to,  36®C.  It  was  determined  in  the 
analysis  that  the  uncertainty  in  these  values  is  on  the  order  of  20  0. 

The  practice  of  relating  all  Ry^,  values  to  the  same  T,,  value  of  36®C  is 
not  necessary  if,  in  the  ATRES/DAAS  software,  all  thermistor  sensors  are 
allowed  to  retain  their  unique  constants  in  the  linear  approximation,  or  if  a 
polynomial  expression  is  used  to  relate  temperature  to  resistance.  In  the 
linear  approximation,  greater  accuracy  is  obtained  when  Ry^  relates  to  the 
midpoint  temperature.  The  common  reference-temperature  is  used  here  in  order 
to  observe  statistical  variations  in  thermistor  constants  as  a  check  on  the 
manufacturer's  stated  tolerances. 

The  determination  of  the  slope,  P,  requires  two  data  points.  Four 
values  of  P  are  determined  in  the  calibration  process;  first  using  the  low 
temperature  T^  and  the  midpoint  temperature  T^,  from  which  we  obtain  P^^,  then 
using  the  midpoint  temperature  and  high  temperature  T^j  giving  a  slope  value 
P,^^  between  these  two  points.  Another  P  value  (P^^)  is  obtained  using  the  low 
and  high  temperature  end  points,  and  finally  an  average  slope  (Pgyg)  is 
calculated  using  the  three  values. 

Each  of  the  p  values  (except  P^^^)  has  a  corresponding  uncertainty  which 
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is  monitored  throughout  the  calibration  process.  If  it  is  assumed  that  a  plot 
of  the  natural  log  of  thermistor  resistance  versus  the  reciprocal  of  absolute 
temperature  is  a  straight  line  over  the  22 °C  to  50 °C  temperature  range,  then 
variations  in  P  values  should  be  within,  or  close  to,  the  predicted 
uncertainties.  Variations  which  fall  outside  predicted  P  uncertainties  can  be 
taken  as  an  indication  either  of  unreliable  temperature/resistance  data  or 
that  the  thermistor  curve  cannot  be  considered  linear  over  the  designated 
temperature  range.  The  maximum  uncertainty  in  p  occurs  for  Pj^j^  and  is 
generally  in  the  vicinity  of  23  Kelvin.  The  most  accurate  p  value  is  Pj^j^  with 
predicted  uncertainties  of  about  8  Kelvin. 

The  development  of  the  calibration  procedure  can  be  summed  up  in  a 
series  of  plots,  which  are  generated  using  an  altered  version  of  program 
ERROR. BAS.  The  plots  show  the  uncertainty  in  T  (thermistor  temperature) 
resulting  from  variations  of  the  measurement  tolerance  (%)  for  P,  T^,  R,  e^, 
and  e^(T).  For  each  plot  the  following  nominal  values  were  used: 

P  =  3500  K 
Wp  =  20  K, 

R^^  =  7,000  n, 

Wrto  =  25  £2, 

Tj,  =  36°C, 

=  0.1‘’C, 
e^  =  2.50  V, 

Wes  =  Weodi)  =  =  0.0001  V. 

The  tolerances  listed  above  were  held  constant  except  where  they  were 
varied  for  a  specific  plot*  The  initial,  later,  and  Delta  temperature  curves 
show  temperature  uncertainty  at  25°C,  50®C,  and  the  difference  (AT)  between 
the  two,  respectively.  In  each  plot  (Figs  4-8),  the  percent  tolerance  on  one 
independent  variable  is  allowed  to  vary  from  zero  to  the  initial  or 
precalibration  tolerance  (for  example,  thermistors  were  purchased  with  a 
manufacturer's  stated  p  tolerance  of  5%  but  are  now  being  calibrated  to  a 
tolerance  of  0.5%)  while  all  other  independent  variables  are  held  constant  at 
the  current  or  calibrated  tolerance.  This  is  summarized  in  Table  11. 
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Variable 


Initial  or 


Current  or 


Plot 


1.2 


Fig.  7.  Temperature  uncertainty  due  to  source- 
voltage  (e^)  tolerance. 


Fig.  8.  Temperature  uncertainty  due  to  voltage 
measurement  (e^  and  e^(T))  tolerance. 

Figure  6  underscores  the  importance  of  having  a  precisely  known,  stable 
circuit  impedance.  Figure  7  was  included  to  show  error  incurred  when  the 
source-voltage  drifts  and  is  not  precisely  monitored  (25%  variations  have  been 
noted  in  the  calibration  process).  This  plot  differs  from  the  next  in  that 


Fig.  8  shows  uncertainty  due  to  the  accuracy  of  all  voltage  measurements 
combined. 

Note  that  no  plot  is  given  for  temperature  uncertainty  versus  thermistor 
resistance  R^.  This  is  because  we  are  not  measuring  resistance  directly,  but 
rather  the  voltage  across  the  voltage  divider  resistor  R.  Any  plot  which 
varies  a  parameter  which  affects  the  accuracy  of  the  R^  values  (e^, 
eg(T„)  and  R),  shows,  indirectly,  the  effect  on  R^. 

Additional  Considerations 

In  addition  to  the  error  sources  and  corrective  terms  already  presented, 
the  effects  of  thermistor  aging  and  instability,  lead-wire  resistance,  direct- 
current  error,  and  system  noise  should  also  be  considered. 

Thermistor  aging  is  a  well  documented  phenomenon^^*^^*^^.  While  it  is 
known  that  the  problem  is  related  to  the  bulk  semiconductor  material  (nickel- 
manganese  oxides),  the  exact  mechanisms  governing  thermistor  instability  due 
to  aging  are  unknown.  Studies  show  a  substantial  increase  in  the  stability  of 
bead-  and  disk-type  thermistors  when  they  are  encapsulated  in  glass  but  the 
reasons  for  this,  too,  are  not  fully  understood. 

A  plot  of  temperature  error  due  to  aging  versus  log  time  will  generally 
be  a  straight  line  for  thermistors  stored  at  constant  temperature.  That  is, 
the  same  change  in  thermistor  characteristics  will  be  encountered  in  the 
period  of  0  to  10  hours  as  in  10  to  100  hours.  Table  12  shows  the  results  of 
an  experiment  in  which  groups  of  5  to  10  thermistors  were  aged— each  group  at 
a  different  constant  temperature — for  periods  of  25  and  5000  hours^^.  Column 
A  gives  the  mean  temperature  shift  for  each  group.  These  should  be  compared 
to  the  predicted  shifts  in  column  B  obtained  by  curve-fitting  the  experimental 
data  to  an  equation  of  the  form 

y=x:*  b. 

where  y  is  the  temperature-shift  in  mK,  x  is  the  hours  of  aging,  and  the 
exponent  a  and  factor  b  are  functions  of  Kelvin  temperature. 
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Table  12 >  Actual  vs.  predicted  shifts  due  to 


Age 

Temoerature 

Time 

Hours 

-80 

25 

-80 

5000 

25 

25 

25 

5000 

70 

25 

70 

5000 

100 

25 

100 

5000 

150 

25 

150 

5000 

200 

25 

200 

5000 

250 

25 

250 

5000 

thermistor-aging 

A 

Actual 

Shift 

mK 

0 

+2 

+1 

+2 

+1 

0 

-1 

-8 

-7 

-25 

-21 

-122 

-32 

-252 


B 

C 

Predicted 

Error 

Shift 

A-B 

mK 

mK 

+8 

-8 

+6 

-4 

+9 

-3 

+7 

-5 

+1 

0 

+3 

-3 

-2 

+1 

-5 

-3 

-9 

+2 

-34 

+9 

-19 

-2 

-106 

-12 

-34 

+2 

-268 

+16 

The  first  several  decades  of  temperature  error  can  be  eliminated  by 
elevating  the  thermistor  to  its  highest  operating  temperature.  This  condition 
is  realized  when  the  thermistors  are  molded  into  the  sensor  disk.  In  this 
step,  the  thermistor  is  brought  from  about  100°C  to  a  maximum  temperature  of 
150 and  back  down  in  a  process  which  takes  about  1  hour. 

Though  the  effects  of  aging  instability  as  presented  in  Table  12  are 
clearly  within  our  predicted  temperature  uncertainty,  they  can  be  conveniently 
monitored  periodically  at  the  triple-point  of  water.  Because,  with  the  use  of 
an  ice  bath,  thermometer  measurements  are  not  necessary,  a  higher  degree  of 
measurement  precision  can  be  obtained,  thereby  bringing  the  variations  into 
our  measurement  realm.  It  is  expected  that  changes  in  the  output  versus 
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temperature  relationship  will  be  apparent  even  outside  the  calibrated 
range^®'^^. 

The  relatively  high  thermistor  impedance  renders  leadwire  resistance 
negligible^®.  For  our  system,  leadwire  resistance  has  been  experimentally 
determined  to  be  in  the  vicinity  of  1  fi — within  the  3  Q  uncertainty  and  hence 
is  not  being  compensated  for. 

When  direct  current  is  run  through  a  thermistor,  errors  can  occur  due  to 
ionic  conduction  in  the  thermistor  and  EMFs  created  at  the  junction  of 
dissimilar  metals  in  the  circuit.  Hajiev  and  Agarunov^^  describe  polarization 
and  electrolysis  processes  which  can  occur  due  to  ionic  conduction  in 
thermistors.  This  can  result  in  irreversible  changes  in  thermistor 
characteristics.  They  report  some  thermistors  which  are  initially  unstable 
and  others  which  become  unstable  after  being  subjected  to  direct  current  for  a 
few  days  or  in  some  instances  a  few  hours.  The  effects  of  thermistor 
instability  due  to  direct  current  can  be  eliminated  by  use  of  an  alternating 
current  bridge  or  by  periodically  reversing  the  current  through  the 
thermistor.  Similarly,  errors  caused  by  thermal  EMFs  within  a  circuit  can  be 
minimized  using  the  same  techniques22,23^ 

While  instability  errors  due  to  direct-current  ionic  conduction  have  not 
been  quantified,  those  due  to  thermal  EMFs  are  expected  to  be  relatively  small 
(for  example  the  EMF  generated  at  a  platinum- copper  junction  is  typically  6-8 
^V/®C)  over  our  temperature  range  of  interest,  and  within  our  predicted 
uncertainty.  However,  in  cases  where  high  precision  is  necessary  the  effects 
of  direct— current  errors  should  be  investigated. 

Because  of  the  high  output  voltage  associated  with  thermometry  (1-2  V 
over  20®C  -  50°C  temperature  range  with  2.5  V  source  voltage),  the  effects  of 
noise  (electrical  and  quantization)  are  expected  to  be  small.  A  study  was 
conducted  where  the  voltage  noise  and  signal-to-noise  ratio  were  determined 
for  three  thermistor/sensors  at  temperatures  of  23®C  and  50°C.  It  was  found 
that  the  average  standard  deviation  of  the  noise  in  a  constant  voltage 
(temperature)  signal  was  0.0018  V  on  a  1.124  V  signal  at  23®C  and  0.0017  V  on 
a  1.78  V  signal  at  50®C.  These  relate  to  signal/noise  ratios  of  600  and  1000, 
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respectively.  After  processing,  the  noise  is  in  terms  of  temperature,  and  is 
experimentally  verified  to  be  in  the  vicinity  of  O.OSOC  at  room  temperature 
for  a  signal/noise  ratio  of  about  500.  Steps  can  be  taken  to  further  reduce 
electrical  noise  by  shielding  cables,  etc.,  and  quantization  noise  can  be 
reduced  by  signal  amplification  prior  to  A/D  conversion,  though  it  is  not 
necessary  at  this  time. 

EXPERIMENTAL  METHOD 
Instrumentation 

Experiments  were  conducted  at  the  PEBr,  FPSD,  U.S.  Army  Natick 
Research,  Development,  and  Engineering  Center,  Natick,  MA.  Figure  9  shows  a 
sketch  of  the  experimental  setup,  followed  by  a  list  and  brief  description  of 
the  equipment  currently  being  used  in  the  thermistor  calibration  procedure. 


Fig.  9.  Experimental  setup  for  thermistor  calibration. 
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Thermometer-1  -  Mercury  total  immersion,  nitrogen  filled.  Range  from 
-10°C  to  lOl^C  in  declinations  of  0.1®C.  Built  by  Arthur  H.  Thomas  Co., 

Phil.,  PA.  Calibrated  by  National  Bureau  of  Standards.  Thermometer-1  is  used 
to  monitor  the  temperature  of  the  bath. 

Thermometer-2  -  Mercury  total  immersion.  Range  from  -0.6°C  to  100. 4°C 
in  declinations  of  0.2°C,  used  to  determine  temperature  of  stem  of 
thermometer-1  for  emergent-stem  correction. 

Power  supply  -  Lambda  Electronics  Corp.,  Melville,  L.I.,  New  York. 

Range  of  0-20  VDC,  used  as  thermistor  power  supply. 

Multimeter  -  Fluke,  Model  8024A  Multimeter,  John  Fluke  Mfg.  Co.,  Inc., 
Seattle,  Washington.  Used  to  set  voltage  supply  with  +0.02  V  accuracy. 

Heating  Mantle  -  Glas-Col  Apparatus  Co.,  Terre  Haute,  Indiana,  Cat.  No. 
TM-108.  Hemispherical  well  filled  with  200-mesh  copper  powder. 

Aluminum  Holder  -  Built  at  Natick,  disk  diameter=3 . 125  inches,  1  inch 
thick.  Drilled  to  hold  four  thermistor/sensor  assemblies  and  mercury  well  of 
Thermometer- 1 . 

Wood  Cover  -  Built  at  Natick,  disk  diameter=6  inches,  1  inch  thick. 
Drilled  to  allow  access  to  thermistor/ sensor  assembly  terminals.  Mounts  on 
top  of  heating  mantle  to  reduce  convective  losses. 

Variac  Autotransformer  -  General  Radio  Company,  Concord,  MA,  type  W5MT3, 
load  0-140  V.  Used  to  control  power  dissipation  in  heating  mantle. 

Switch  -  12  position,  used  to  select  thermistor/ sensor  assembly  under 

study. 
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Linear  Voltage  Divider  -  Built  at  Natick.  Contains  single  10  K  resistor 
calibrated  at  Natick.  True  resistance  at  24.8  C  =  9998.8  0  ±0.002%,  0±50 
ppm/^C. 


Hardware  - 

1.  Zenith  Data  Systems  Model  ZWX-0248-62  PC  -  Heath 
Company,  Benton  Harbor,  Michigan.  Used  in  data  acquisition  and 
processing. 

2.  Dash-8  -  Data  Acquisition  and  Control  Interface  Board 
installed  in  Zenith  PC  to  provide  analog/digital  interface  and 
high-speed  data  acquisition. 

3.  EXP-16  -  Expansion  Multiplexer  and  Instrumentation 
Amplifier  used  to  extend  Dash-8  capabilities.  Both  Dash-8  and 
EXP-16  are  products  of  MetraByte  Corporation,  Taunton,  MA. 

Software  - 

1.  Labtech  Notebook^”  -  Software  package  from  Laboratory 
Technologies  Corporation,  Wilmington,  MA.  Used  in  data 
acquisition.  Interfaces  with  Dash-8  and  EXP-16. 

2.  Lotus  1-2-3™  -  LOTUS,  Cambridge,  MA.  Used  to 
process  data  files  generated  with  Labtech  Notebook. 

Experimental  Procedure 


1.  Labtech  Notebook  is  set  up  for  two  input  channels;  one  to  monitor 
source  voltage  and  the  other  for  voltage-divider  output -volt age.  Sampling 
rate  is  set  at  50  Hz  (normal  speed)  for  a  duration  of  10  s.  Both  channels  are 
dumped  into  one  file,  which  is  placed  in  the  1—2—3  directory. 

2.  Four  thermistor/ sensor  assemblies  are  placed  in  the  aluminum  holder 
which  is  buried  in  the  200-mesh  copper  powder  inside  the  well  of  the  heating 
mantle.  A  cardboard  cover  is  placed  over  the  aluminum  holder  as  an  additional 
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insulator  and  the  wooden  cover  is  placed  over  the  whole  assembly  such  that  the 
sensor  terminals  can  be  accessed  through  the  drilled  holes. 


3.  Four  pairs  of  alligator  clips  are  attached  to  the  terminals  taking 
special  note  that  the  clip  numbers  correspond  to  the  proper  location  of  the 
sensor  in  the  A1  holder.  Plastic  funnels  are  then  fitted  over  the  alligator 
clips  and  the  foam  plugs  are  inserted  into  the  funnel  openings  to  reduce 
convective  losses, 

4.  The  12-position  switch  is  set  to  a  neutral  position  and  the  source- 
voltage  is  turned  on  and  adjusted  with  the  Multimeter  to  an  output  of  0.50  V. 
(Note:  Allowing  a  current  to  run  through  a  sensor  will  cause  it  to  heat  up  and 
prevent  the  assembly  from  reaching  thermal  equilibrium. ) 

5.  An  O-ring  is  slid  over  the  mercury  well  of  thermometer-1  and 
positioned  where  the  glass  meets  the  mercury  well.  The  O-ring  should  be 
positioned  such  that  it  rests  on  top  of  the  wood  insulator  plate.  This 
reduces  convective  losses  through  the  accommodating  hole  and  also  prevents  the 
mercury  well  from  resting  on  the  bottom  of  the  Al  holder.  (Note:  Allowing  the 
thermometer  bulb  to  contact  the  Al  holder  in  this  vertical  position  will 
create  pressure  in  the  bulb  and  a  bias  in  the  temperature  reading.) 

6.  Thermometer-2  is  attached  to  thermometer-1  such  that  the  two  are  in 
contact  along  the  length  of  the  stem.  The  system  is  now  allowed  to  sit  for  20 
-  30  min  to  allow  the  assembly  to  reach  thermal  equilibrium.  It  is  important 
that  the  room  temperature  remains  stable  within  this  period  because  the  stem 
of  thermometer-1  is  acting  as  a  heat-transfer  fin.  Once  the  assembly  reaches 
thermal  equilibrium  (evidenced  by  no  change  in  temperature  indicated  by 
thermometer-1),  the  low-power  room-temperature  data  can  be  taken. 

7.  The  switch  is  positioned  such  that  current  is  running  through  the 
sensor  in  location-1  of  the  Al  holder.  Acquisition  of  data  is  begun  with 
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Labtech  Notebook.  The  switch  is  left  in  position-1  for  two  seconds,  then  rn 
position-2  for  two  seconds,  etc.,  until  data  have  been  collected  for  all  four 
sensors,  after  which  the  switch  should  again  be  set  to  a  neutral  positron. 

8.  Labtech  Notebook  is  exited  and  Lotus  1-2-3  is  entered,  retrieving 
file  CALl.WKl.  The  Labtech  data  file  is  imported  into  the  appropriate 
spreadsheet  column  and  the  file  is  saved  under  a  new  name,  making  note  of  the 
calibration  batch  (for  example,  BATCHl.WKl). 

9.  Returning  to  Labtech  Notebook,  the  source-voltage  is  changed  to  2.50 
V  using  the  Multimeter  to  monitor  the  output. 

10.  It  is  important  to  observe  whether  the  temperature  has  changed 
according  to  thermometer-1.  (A  0.05»C  change  is  acceptable,  but  a  change  of 
a.l^C  or  greater  is  not,  and  the  complete  process  from  Step  6  must  be 
repeated.)  Step  7  must  be  repeated  at  the  new  voltage  setting  and  the  Labtech 
Notebook  data  must  be  entered  into  the  appropriate  column  of  the  new  Lotus 
file. 


11.  The  temperature  should  be  recorded  in  the  Lotus  spreadsheet  and 
also  on  the  ralibration  Sheet.  Lotus  will  automatically  average  the  low-  and 
high-power  source-voltage  (columns  A  and  E).  These  must  be  entered  into  the 
Calibration  Sheet.  The  voltage  corresponding  to  each  sensor  must  be  averaged 
using  Lotus,  and  these  values  entered  in  the  Calibration  Sheet  for  both  low- 
and  high-power. 

12.  The  variac  is  turned  on  at  a  setting  of  30.  The  heating  mantle 
should  be  allowed  to  increase  in  temperature  until  thermometer-1  reads  28®C. 
At  this  point  the  variac  is  turned  off  and  the  assembly  is  allowed  to  reach 
thermal  equilibrium.  The  system  should  equilibrate  close  to  36°C — the 
midpoint  temperature. 
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13. 


Thermometer-1  must  be  observed  closely,  and  the  voltage  data 
recorded  as  described  above  at  the  point  where  the  assembly  just  reaches 
thermal  equilibrium  close  to  36®C.  The  stem  temperature  must  also  be  recorded 
as  indicated  by  thermometer-2.  (In  general,  thermometer  readings  are  always 
taken  as  the  temperature  moves  from  low  to  high,  because  this  is  the  way  in 
which  thermometers  are  calibrated  both  by  the  manufacturers  and  by  the 
National  Bureau  of  Standards.) 

14.  After  completion  of  Step  13,  the  variac  is  again  turned  on  and  left 
on  until  thermometer-1  reads  43.5»C.  The  system  will  reach  equilibrium  close 
to  50 “C.  The  necessary  data  are  recorded  as  with  the  midpoint  temperature. 
There  is  now  the  necessary  information  to  calculate  thermistor  constants 
P ,  and  5 . 


15,  Thermistor  constants  are  calculated  using  the  QuickBASIC^^  program 
entitled  DC_BETA,BAS.  When  run,  the  program  prompts  for  the  following 
information: 

TL  -  Low  temperature  obtained  with  thermometer-1. 

VSL  -  Low-temperature  high-power  source-voltage. 

LOW  POWER  VSL  -  Low- temperature  low-power  source-voltage. 

TM  -  Midpoint-temperature  obtained  with  thermometer-1. 

VSM  -  Midpoint-temperature  source-voltage. 

STEM  TEMP  (M)  -  Stem-temperature  at  midpoint. 

TH  -  High-temperature  obtained  with  thermometer-1 
VSH  -  High-temperature  source-voltage. 

STEM  TEMP  (H)  -  Stem-temperature  at  high  point. 

VTL  -  Voltage-divider  voltage  at  low-temperature  high-power. 

LOW  POWER  VTL  -  Voltage-divider  voltage  at  low-temperature  low-power. 

VTM  -  Voltage-divider  voltage  at  midpoint -temp. 
yipH  -  Voltage-divider  voltage  at  high-temp. 

The  system  iterates  on  the  dissipation  constant  then  gives  the  output  as 
in  the  following  example: 
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Table  13.  Output?  of  program  DC  BETA. BAS _ for — sensor  T97. 


BETA,  Kelvin 

Uncertainty 

RT,  ohms 

BLM 

3496.13 

20.08985 

12525.42 

BMH 

3497.361 

18.16944 

7790.396 

BLH 

3496.759 

9.629001 

4742.007 

BAVE 

3496.75 

TEMP 

EXP  TEMP 

ERROR 

LOW 

24.02902 

24.0306 

-1.577377E-03 

MED 

36.59214 

HIGH 

50.76709 

50.7698 

-2.708435E-03 

RT36 

7942.424 

DC 

2.948801E-04 

SHL 

SHM 

SHH 

0.5290134 

0.5266719 

0.4669011 

Given  under  the  heading  'BETA'  are  the  four  p  values  and  in  the  next 

column  labeled  'UNCERTAINTY'  is  the  uncertainty  in  units  of  Kelvin  for  each  of 

the  respective  p  values.  These  should  be  compared  to  the  'BAVE'  value.  The 

three  resistance  values  listed  under  'RT'  correspond  to  the  three  temperatures 

listed  under  'TEMP'.  The  temperatures  listed  under  'TEMP'  are  the  'actual' 

thermistor  temperatures  obtained  by  adding  the  self-heating  values  listed 

under  'SHL',  'SHM',  and  'SHH',  and  the  emergent-stem  corrections  to  the 

temperatures  obtained  with  thermometer-1  during  the  calibration  process.  The 

temperatures  given  under  'EXP  TEMP'  are  the  temperatures  calculated  using  the 

linear  approximation,  the  listed  RT  values,  and  the  experimentally  determined 

B  and  R,  values.  The  'ERROR'  values  are  obtained  by  subtracting  the  'EXP 
^^ave  To 

TEMP'  listings  from  the  'TEMP'  listings  and  are  included  as  a  check  they  have 
no  statistical  significance.  The  thermistor  reference-temperature  resistance 
R  relating  to  T  =36'»C  is  given  as  'RT36',  and  the  dissipation  constant  is 
given  as  'DC'*  All  of  these  values  should  be  entered  into  the  Thermistor  Data 
Sheet  *  This  completes  the  thermistor  calibration  procedure. 
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RESULTS 


The  uncertainty  analysis  predicts  a  standard  deviation  between  O-l^C  and 
0.2®C  over  the  temperature  range  of  22 °C  to  50“C.  Results  to  date  show 
absolute  variations  of  0.1 °C  or  less.  At  the  writing  of  this  report,  64 
sensors  have  been  calibrated  as  described  above.  In  order  to  quantify  the 
precision  with  which  the  thermistors  have  been  calibrated,  three  different 
studies  were  performed.  In  each  study,  three  or  four  sensors  are  compared. 

The  inputs  range  from  20»C  to  SO’C  in  steps  of  a.l^C.  Thermistor 
resistance  is  calculated  over  this  temperature  range  for  the  group  using  each 
sensor's  unique  p,  and  5  value.  The  dissipation  constant  is  implemented 

such  that  thermistor  self-heating  is  subtracted  from  the  temperatures  from 
which  the  Ry  values  were  calculated  (2.50  V  source-voltage  is  assumed).  The 
three  or  four  temperature  curves  are  then  averaged  and  each  curve  is 
subtracted  from  the  average  creating  absolute-error  curves  which  are  plotted 
on  a  single  graph  for  comparison. 

In  the  first  study,  six  sensors,  chosen  at  random  from  different 
batches,  are  calibrated,  three  times  each,  to  check  for  consistency  in 
results.  The  curves  from  each  of  the  three  calibration  runs  are  then  plotted. 


Fig.  10.  Absolute-error  curves  for  three  calibration  runs 

for  sensor  T76. 
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The  maximum  deviation  from  the  average 


as  described  above,  for  each  sensor, 
for  any  one  sensor  occured  with  sensor  T76  and  was  about  0.085®C.  This  is 

shown  in  Fig  10. 

The  smallest  deviation  from  the  average  occurred  with  sensor  T80  and  was 
about  0.028®C.  This  is  shown  in  Fig.  11. 


Fig.  11.  Absolute-error  curves  for  three  calibration  runs 

for  sensor  T80. 

For  the  second  study,  19  batches  of  four  sensors  each  were  plotted  as 
described  above  in  order  to  observe  variations  within  a  batch.  For  each  of 
the  19  batches,  absolute  deviations  are  within  0.1°C  from  average.  Figure  12 
shows  the  error  curves  for  Batch-19,  which  exhibited  the  greatest  deviation 

from  average. 
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Fig.  12.  Absolute-error  curves  for  sensors  in  Batch-19. 
An  exceptionally  precise  batch,  Batch-5A,  is  shown  in  Fig.  13 


Fig.  13.  Absolute-error  curves  for  sensors  in  Batch-5A. 

For  the  third  study,  individual  sensors  were  randomly  selected  one  from 
each  batch“and  grouped  into  four  batch-mixes.  This  was  done  in  order  to 
observe  variations  from  one  batch  to  another.  For  each  batch-mix,  absolute 
deviations  are  less  than  0.1°C  and  within  predicted  uncertainty.  An  example 
is  shown  in  Fig*  14. 
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Fig.  14.  Absolute-error  curves  for  Batch-mix  3. 

It  should  be  noted  that  in  the  above  three  studies  it  is  assumed  that 
there  is  a  random  scatter  of  data,  and  that  the  average  of  a  number  of  curves 
is  the  most  accurate.  The  consistency  of  the  results  would  seem  to  support 
this  assumption.  However,  it  is  possible  that  the  entire  data  base  is  biased, 
and  the  random  variations  are  superimposed  over  this  bias.  Because  we  have  no 
additional  standard  against  which  to  verify  our  results  (i.e.,  a  factory- 
calibrated  thermistor),  there  is  no  way  of  knowing.  However,  because  the 
current  standard  is  a  NBS-calibrated  thermometer,  it  is  believed  that  rf  there 
is  such  a  bias,  it  is  probably  quite  small  and  presumably  within  predicted 
uncertainty.  Though  we  cannot  at  the  present  time  verify  high  accuracy,  we 
have  demonstrated  high  precision. 

In  developing  the  software  to  be  used  in  data  acquisition  and  processing, 
it  is  necessary  to  know  the  magnitude  of  variations  in  thermistor  constants. 
Recall  that  the  thermistors  were  purchased  with  an  tolerance  of  ±25% 
relative  to  10,000  £2  at  T^j=25‘‘C.  Also,  P  values  were  expected  to  vary,  for  a 
given  batch,  by  ±3%  from  some  nominal  value,  and  for  different  batches,  by  as 
much  as  5%.  From  the  64  thermistor/ sensors  calibrated  to  date,  the  following 

information  was  obtained; 

Reference-Temperature  Resistance  R^^  (36®C); 

Average  =  6736  £2 

Standard  Deviation  =  564  £2  (8%) 
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Maximum  Value  =  10213  Q 

Minimum  Value  =  6193  Q 

%  Maximum  Deviation  >  Average  =  52% 

%  Maximum  Deviation  <  Average  =  8% 

The  large  %  maximum  deviation  for  values  greater  than  average  is 
evidence  of  the  existence  of  an  outlier,  verified  to  be  the  maximum  value 
shown  above*  This  unusually  high  Rj^  value  was  double-checked  by  calibrating 
the  thermistor  a  second  time,  and  was  found  to  be  correct.  When  the  outlier 
is  removed,  however,  the  variations  fall  within  the  manufacturer's  stated 
tolerance  of  ±  25%  as  shown  below: 

Average  =  6681 

Standard  Deviation  =  355  (5%) 

Maximum  Value  =  7942  Q 
Minimum  Value  =  6193  Q 
%  Maximum  Deviation  >  Average  =  18% 

%  Maximum  Deviation  <  Average  =  8% 

A  similar  study  was  conducted  for  variations  in  with  the  following 

results: 

Beta: 

Average  =  3467  Kelvin 
Standard  Deviation  =  96  K  (2.8%) 

Maximum  Value  =  3620  K 
Minimum  Value  =  3119  K 
%  Maximum  Deviation  >  Average  =4.4% 

%  Maximum  Deviation  <  Average  =  10% 

The  minimum  value  is  an  outlier,  when  it  is  removed  we  obtain: 

Average  =  3472  Kelvin 
Standard  Deviation  =  85  K  (2.4%) 

Maximum  Value  =  3620  K 
Minimum  Value  =  3251  K 
%  Maximum  Deviation  >  Average  =  4.4% 

%  Maximum  Deviation  <  Average  =  6.2% 
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The  deviations — both  greater  and  lesser  than  average — are  outside  the 
manufacturer's  stated  3%  tolerance.  Tolerance  variations  for  both  and  p 
have  little  consequence  if  each  sensor  is  allowed  to  retain  its  own 
unique  constants  in  the  linear  approximation,  other  than  the  fact  that  P  will 
be  linear  over  a  different  temperature  range. 

In  the  experimental  determination  of  the  dissipation  constant  it  has 
been  noticed  that  approximately  1  of  8  thermistor/sensors  gives  a  dissipation 
constant,  which  is  unusually  lower  or  higher  than  expected.  For  the  64 
sensors  calibrated  to  date,  the  following  information  was  obtained: 

Average  5  =  2.868x10'^  W/^C 
Standard  Deviation  =  0.534x10'^  W/®C 
Predicted  Uncertainty  =  0.23x10'^  W/®C 
Though  it  is  expected  that  there  will  be  variations  in  5  due  to  varying 
heat-transfer  mechanisms,  these  larger-than-expected  fluctuations  are  altering 
the  self-heating  corrective  term  such  that  the  corrected  temperatures  are  as 
much  as  0.2  -  0.3®C  in  error  with  the  rest  of  the  sensors  in  the  batch.  When 
this  occurs,  the  sensor  is  often  recalibrated,  usually  with  the  same  result. 
The  exact  cause  of  this  is  not  known,  but  a  nonlinear  power/ conductance 
relationship  in  the  conductive  epoxy  used  to  electrically  connect  the 
thermistors  leads  to  the  terminals  has  been  considered.  No  literature  has 
been  found  to  support  this  hypothesis,  however.  Thermistor  instability  is 
another  possibility,  but  neither  has  this  been  confirmed. 

The  current  method  of  dealing  with  this  discrepancy  is  by  implementing 
an  average  dissipation-constant  calculated  from  existing  data  and  given  above 
as  2.867x10'^  W/®C.  Program  DC_BETA.BAS  was  altered  to  use  this  average  5 
value.  The  new  program,  AVG_DCB.BAS,  does  not  iterate  on  5  and  p,  but  rather 
calculates  thermistor  constants  using  the  average  6  value.  Figure  15  shows 
absolute-error  curves  for  Batch-19  for  which  the  average  5  value  was 
implemented.  Note  that  the  error  curves  are  all  tied  down  at  eibout  24®C,  the 
temperature  at  which  the  individual  unique  dissipation  constants  were 
calculated  for  this  batch. 
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Fig.  15.  Absolute-error  curves  for  Batch-11  for 
which  an  average  5  was  used. 

The  high  precision  implied  by  this  error  plot  is  supportive  of  the  use 
of  an  average  6.  Additional  studies  may  have  to  be  taken  along  these  lines  in 
order  to  prove  the  validity  of  this  approach. 

CONCLUSIONS  AND  RECOMMENDATIONS 

A  procedure  has  been  developed  and  utilized  for  the  calibration  of 
thermistors  mounted  in  simulated  skin  sensors.  An  uncertainty  analysis  has 
been  performed  based  on  a  linear  approximation  of  the  temperature/resistance 
characteristics  of  negative— temperature— coefficient  thermistors.  Expressions 
were  derived,  and  uncertainties  predicted  for  both  absolute-  and  relative- 
temperature  measurements. 

In  the  analysis,  a  number  of  parametric  measurement  schemes  were 
considered.  The  accuracy  of  temperature  measurements  made  with  a  NBS 
calibrated  thermometer  was  compared  to  those  made  with  a  precalibrated 
thermistor  with  the  result  that  thermometer  measurements  are  found  to  be  the 
more  accurate.  Resistance  values  obtained  through  voltage  measurements  were 
found  to  have  a  lesser  degree  of  uncertainty  than  those  made  with  an  available 
ohmmeter.  Through  the  experimental  determination  of  P  and  the  analysis 
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predicts  a  temperature  uncertainty  of  0.1-0. 2 °C  over  the  temperature  range  of 
25  -  50®C. 

In  addition,  a  measurement  scheme  was  devised  for  the  experimental 
determination  of  thermistor  self-heating  due  to  power  dissipation,  and 
temperature  correction  terms  have  been  implemented  for  this  and  thermometer 
emergent-stem  errors. 

Software  was  written  which  calculates  thermistor  constants  using  voltage 
and  temperature  measurements  while  simultaneously  predicting  uncertainties  in 
the  calculated  values,  so  that  the  integrity  of  each  calibration  run  can  be 
closely  monitored. 

Results  from  the  64  sensors  calibrated  as  thermistors  to  date  show 
variations  of  less  than  a.l^C;  within  the  predicted  temperature-measurement 
uncertainty  and  required  ±0.5®C  accuracy. 

A  decision  will  have  to  be  made  concerning  the  form  of  the  expression 
used  to  relate  temperature  to  resistance  in  the  ATRES/DAAS  data  processing 
software.  The  most  precise  temperature/resistance  relationship  is  obtained 
through  the  polynomial  expression  (Eq.  7).  The  constants  a^,  a^,  and  a^  are 
not  currently  being  determined  for  each  thermistor/ sensor.  The  software  is, 
however,  setup  for  this  purpose  and  it  would  not  require  much  time  to  complete 
the  task.  In  this  case,  the  self-heating  correction  term  must  still  be 
utilized  in  order  to  account  for  variations  in  source  voltage,  and  hence, 
variations  in  dissipated  power. 

The  error  (±0.1®C)  inherent  in  the  linear  approximation  due  to  nonlinear 
P  is  almost  negligible  when  one  considers  that  our  temperature-measurement 
uncertainty  is  between  0.1  and  0.2®C.  For  this  reason,  it  is  believed  that 
the  linear  approximation  represents  an  adequate  relationship  between 
temperature  and  resistance  for  our  purposes. 

If  processing  speed  is  a  consideration  (real  time  display),  the  linear 
approximation  can  be  used  with  average  p,  and  6. 

In  the  event  that  greater  precision  is  required  for  a  particular 
application,  a  higher  degree  of  accuracy  can  be  obtained  with  a  more  precise 
ohmmeter  for  resistance  measurements  and  a  factory-calibrated  thermistor  for 
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temperature  measurements.  Essentially  all  of  the  software  developed  for  this 
calibration  procedure  is  currently  applicable,  or  can  be  altered  to 
accommodate  a  more  precise  calibration  method. 

It  is  believed  that  the  calibration  process  presented  here  represents  a 
reliable  technique  by  which  to  achieve  the  required  temperature-measurement 
precision  of  the  new  ATRES/DAAS  system  currently  being  developed  at  Natick. 
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APPENDIX  COMPUTER  PROGRAMS 


!•  Program  ERROR. BAS 

2 .  Program  RTNOT . BAS 

3 .  Program  RTNOTA . BAS 

4.  Program  DC_BETA.BAS 
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1.  Program  ERROR. BAS 


5  REM  Program  ERROR. BAS  Brian  Kimball,  2/19/92  ^  ,  ,4.  n, 

10  REM  Calculates  uncertainty  in  absolute  and  delta  i 
measurements 

20  INPUT  "BETA,  KELVIN  B=",  B 

30  INPUT  "Reference  temperature,  celsius  Tr— ",  Tr 
40  INPUT  "Initial  temperature,  celsius  Ti=",  Ti 
50  INPUT  "Later  temperature,  celsius  Tn=",  Tn 
60  REM  INPUT  "Resistor,  ohms  R=",  R 
62  R  =  9998.8  ^ 

65  INPUT  "Thermistor  resistance  at  reference  temperature,  ohms 


70 

80 

90 


RTr=",  RTr 

INPUT  "Source  voltage,  volts  es-“,  es 
INPUT  "Uncertainty  in  Beta,  Kelvin  WB=",  WB 

INPUT  "Uncertainty  in  reference  temperature  resistance,  ohms 


WRTr 

100  REM  INPUT  "Uncertainty  in  resistor  resistance,  ohms  WR=",  WR 

120  INPUT  "Uncertainty  in  reference  temperature,  celsius  WTr-", 
WTr 


125  REM 

130  Tr  =  Tr  +  273.15 

135  Ti  =  Ti  +  273.15 

140  Tn  =  Tn  +  273.15  ,  m 

150  RTi  =  RTr  *  EXP(B  *  ((1  /  Tl)  -  (1  /  Tr)  )  ) 

155  Wes  =  es  *  .0001 
160  eri  =  R  *  es  /  (RTi  +  R) 

170  Weri  =  eri  *  .0001 

180  RTn  =  RTr  *  EXP(B  *  ((1  /  Tn)  -  (1  /  Tr) ) ) 

190  ern  =  R  *  es  /  (RTn  +  R) 

200  Wern  =  ern  *  .0001 
205  REM 

210  REM  PARTI ALS  TO  FIND  ERROR  IN  RTl  AND  RTn 
220  Pli  =  ((es  /  eri)  -  1) 

230  P2i  =  R  /  eri 

240  P3i  =  R  *  es  /  (eri  2)  .  .  ^  ^ 

250  WRTi  =  ( (Pli  *  WR)  ''  2  +  (P2i  *  Wes)  2  + 
^  .5 

255  REM  PRINT  "WRTi=" ,  WRTi 
260  Pin  =  ((es  /  ern)  -  1) 

270  P2n  =  R  /  ern 

280  P3n  =  R  *  es  /  (ern  2) 

290  WRTn  =  ((Pin  *  WR)  "  2  +  (P2n  *  Wes)  ^  2  + 
"  .5 


(P3i  *  Weri) 


(P3n  *  Wern) 


-  2) 


-  2) 


295  REM 

299  REM 

300  REM  PARTIALS  TO  FIND  ABSOLUTE  ERROR  IN  Ti  AND  Tn 
310  P4i  =  (B  *  Tr  "  2)  /  (RTi  *  (Tr  *  LOG (RTi  /  RTr) 


+  B) 


-  2) 


2)  *  LOG(RTi  /  RTr)  /  ( (Tr  *  LOG(RTi  /  RTr)  +  B) 


330  P6i  =  (Tr 
2 ) 

P7i  =  B  "  2  /  ((Tr  *  LOG(RTi  /  RTr)  +  B)  ^  2) 

340  WTi  =  ((P4i  *  WRTi)  2  +  (P5i  *  WRTr)  "  2  +  (P6l  *  WB) 
(P7i  *  WTr)  ^2)  ^  .5 


2  + 


350  REM 
360  P4n  =  (B  *  Tr 
370  P5n 


2) 

2  +  (P6n  *  WB)  ^  2  + 


ID  "  xa.  2)  /  (RTn  *  (Tr  *  LOG(RTn  /  RTr)  +  B)  ^  2) 

TTn  *  Tr  ^  2)  I  (RTr  *  (Tr  *  LOG(RTn  /  RTr)  +  B)  2) 

380  P6n  =  (Tr  "  2)  *  LOG(RTn  /  RTr)  /  ( (Tr  *  LOG(RTn  /  RTr)  +  B) 
^2) 

385  P7n  =  B  -  2  /  ((Tr  *  LOG (RTn  /  RTr)  +  B) 

390  WTn  =  ((P4n  *  WRTn)  2  +  (P5n  *  WRTr) 

(P7n  *  WTr)  ''2)  .5 

PPM  PARTIALS  TO  CALCULATE  UNCERTAINTY  IN  DELTA  T 
420  De55  *  LOG(R?n  /  RTr)  +  B)  *  (Tr  *  LOG(RTi  /  RTr)  +  B )  ^ 

430  P7  =  (((DEN)  *  (Tr  ^  2  *  LOG(RTi  /  RTn)))  ”  ( (B  Tr  2  ^ 
LOG(RTi  /  RTn))  *  (2  *  B  +  Tr  *  (LOG((RTi  *  RTn)  /  (R  *  Tr 

440  P8  =’((('Drar**(2’  *  B  *  LOG(RTi  /  RTn)))  -  ( (B  *  Tr  "  2  * 
MG(RTi  /  RTn))  i  (2  *  Tr  ♦  LOG(RTn  /  RTr)  *  LOG(RTl  /  RTr)  + 
B  *  LOG((RTi  *  RTn)  /  (RTr  ^  2)))))  /  (DEN  2) 

450  P9  =  -((B  *  Tr  "  2  *  LOG(RTi  /  RTn))  *  ( (Tr  /  RTr)  *  (2  *  B  + 
Tr  *  LOG((RTn  *  RTi)  /  (RTr  ^  2)))))  /  (DEN  2)^ 

460  PIO  =  ((DEN)  *  (B  *  (Tr  "  2)  /  4^  .^2^^ 

LOG(RTi  /  RTn))  *  (Tr  *  LOG(RTn  /  RTr)  +  B) )  /  (DEN  2) 

470  Pll  =  ((DEN)  *  (-B  *  (Tr  ^  2)  /  RTn)  -  / 


*  ^  *  MmJ  f  f  f  \  / 

-  ( (B  *  Tr  ^  3  *  LOG (RTi 


/  RTn)  /  RTn)  *  (Tr  *  LOG (RTi  /  RTn)  +  B) ) )  /  (DEN  "  2) 

/  KJ-Ii;  /  .  .  _  .  '  J.  T.Tin^\  X  O  a.  /DO  *  WPTT-^ 


2  +  (P9  *  WRTr)  ^  2  + 

.5 


480  WDelT  =  (  (P7  *  WB)  2  +  (P8  *  WTr) 

(PIO  *  WRTi)  ^  2  +  (Pll  *  WRTn)  "  2) 

485  PRINT 
490  REM 

500  REM  PRINT  RESULTS  ^  .  4.  ••  4-4=1  4.^„,r^  Ti  ic^ii 

510  PRINT  ”The  uncertainty  in  the  absolute  initial  temp  / 

520  PRINT  "The  uncertainty  in  the  absolute  later  temp  Tn  is”,  WTn 

530  PRINT  "The  uncertainty  in  the  Delta-T  measurement  is  ,  W  e 

540  PRINT  "The  uncertainty  in  Delta-T  due  to  the  uncertainty  in:" 

550  WP7  =  P7  *  WB 
560  PRINT  "  Beta",  WP7 
570  WP8  =  P8  *  WTr 
580  PRINT  "  Tr",  WPS 
590  WP9  =  P9  *  WRTr 
600  PRINT  "  RTr",  WP9 
610  WPIO  =  PIO  *  WRTi 
620  PRINT  "  RTi",  WPIO 
630  WPll  =  Pll  *  WRTn 
640  PRINT  "  RTn",  WPll 

650  PRINT  .  4.  •  4.  •  II 

660  PRINT  "Uncertainty  in  Ti  due  to  uncertainty  in: 

670  PRINT  "  Beta",  (P6i  *  WB) 
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680  PRINT 
690  PRINT 
700  PRINT 
710  PRINT 
720  PRINT 
730  PRINT 
740  PRINT 
750  PRINT 
760  PRINT 
770  PRINT 
800  END 


•I  Tr”,  (P7i  *  WTr) 

II  RTr”,  (P5i  *  WRTr) 

”  RTi”,  (P4i  *  WRTi) 

"Uncertainty  in  Tn  due  to  uncertainty  in: 
"  Beta",  (P6n  *  WB) 

II  Tr",  (P7n  *  WTr) 

"  RTr",  (P5n  *  WRTr) 

"  RTn",  (P4n  *  WRTn) 
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2.  Program  RTNOT.BA8 


10  REM  program  RTNOT.BAS,  Brian  Kimball,  2/20/92 

15  REM  This  program  computes  values  of  RTo.  It  also 

20  REM  calculates  uncertainty  in  experimentally  determined 

30  REM  values  of  RTo.  This  version  used  for  calibrating 

31  REM  reference  thermistors. 

40  INPUT  "Ambient  temperature,  Celsius  T=”,  T 

41  IF  T  =  0  THEN 

42  GOTO  210 

43  ELSE 

44  GOTO  50 

45  END  IF 

50  REM  INPUT  "Uncertainty  in  ambient  temperature,  Celsius  WT=", 
WT 

55  WT  =  .05 

60  INPUT  "Thermistor  resistance  at  ambient  temperature,  ohms 
RT=",  RT 

70  REM  INPUT  "Uncertainty  in  RT,  ohms  WRT=",  WRT 

75  WRT  =  3 

80  T  =  T  +  273.15 

90  B  =  3474 

100  WB  =  18 

110  Tr  =  298.15 

111  RTr  =  RT  /  (EXP(B  *  (1  /  T  -  1  /  Tr))) 

120  REM  Calculate  partials 

130  PI  =  1  /  (EXP(B  *  (1  /  T  -  1  /  Tr))) 

140  P2  =  RT  *  (1  /  T  -  1  /  Tr)  /  (EXP(B  *  (1  /  T  -  1  /  Tr))) 

150  P3  =  RT  *  B  /  ((T  "  2)  *  (EXP(B  *  (1  /  T  -  1  /  Tr) ) ) ) 

160  WRTr  =  ((PI  *  WRT)  ^  2  +  (P2  *  WB)  ''  2  +  (P3  *  WT)  ^2)  ''  .5 

165  PEI  =  ((RTr  -  10000)  /  10000)  *  100 
167  PE2  =  (WRTr  /  RTr)  *  100 

169  PRINT 

170  PRINT  "Reference  temp  resistance  is",  RTr 
180  PRINT  "Thermistor  tolerance  is,  percent",  PEI 

190  PRINT  "The  uncertainty  in  this  measurement  (RTr)  is,  ohms", 
WRTr 

200  PRINT  "or,  percent",  PE2 

201  PRINT  "Resistance  term";  (PI  *  WRT) 

202  PRINT  "Beta  term";  (P2  *  WB) 

203  PRINT  "Temperature  term";  (P3  *  WT) 

204  PRINT 

205  PRINT 

206  GOTO  40 
210  END 


3.  Program  RTNOTA.BAS 


10  REM  program  RTNOTA.bas,  Brian  Kimball,  2/24/92 
20  REM  This  program  calculates  the  reference  temperature 

30  REM  resistance  of  a  thermistor.  It  calculates  the  ambient 

31  REM  temperature  by  using  the  resistance  value  of  a 

32  REM  pre-calibrated  reference  thermistor  (T51).It  uses  this 

33  REM  temperature  and  resistance  value  of  the  test  thermistor 

34  REM  (input)  at  this  ambient  temperature  to  determine  the 

35  REM  reference  temperature  resistance  of  the  test  thermistor. 

36  REM  The  reference  temperature  used  here  is  25  degrees  C. 

40  INPUT  "Reference  thermistor  resistance,  ohms  RT51=”,  RT51 

41  IF  RT51  =  0  THEN 

42  GOTO  210 

43  ELSE 

44  GOTO  46 

45  END  IF 

46  INPUT  "Test  thermistor  resistance,  ohms  RT=",  RT 
55  WT  =  .25 

75  WRT  =  20 
80  B  =  3474 
85  WB  =  104 

90  TR  =  298.15 

91  RT51r  =10316 

92  T  =  B  *  TR  /  (TR  *  L0G(RT51  /  RT51r)  +  B) 

115  RTr  =  RT  /  (EXP(B  *  (1  /  T  -  1  /  TR))) 


120  REM  Calculate  partials 

130  PI  =  1  /  (EXP(B  *  (1  /  T  -  1  /  TR))) 

140  P2  =  RT  *  (1  /  T  “  1  /  TR)  /  (EXP(B  *  (1  /  T  -  1  /  TR))) 

150  P3  =  RT  *  B  /  ((T  ^  2)  *  (EXP(B  *  (1  /  T  -  1  /  TR) ) ) ) 

160  WRTr  =  ((PI  *  WRT)  2  +  (P2  *  WB)  ^  2  +  (P3  *  WT)  2)  ^  .5 

165  PEI  =  ( (RTr  -  10000)  /  10000)  *  100 
167  PE2  =  (WRTr  /  RTr)  *  100 


168  PRINT 

169  PRINT 

170  PRINT 
180  PRINT 
190  PRINT 

200  PRINT 

201  PRINT 

202  PRINT 

203  PRINT 

204  PRINT 


"Ambient  temp,  Celsius",  (T  -  273.15) 
"Reference  temp  resistance  is",  RTr 
"Thermistor  tolerance  is,  percent",  PEI 
"The  uncertainty  in  this  measurement  (RTr)  is, 
WRTr 

"or,  percent",  PE2 
"Resistance  term";  (PI  *  WRT) 

"Beta  term";  (P2  *  WB) 

"Temperature  term";  (P3  *  WT) 


ohms" , 
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205  PRINT 

206  GOTO  40 
210  END 


4.  Progrzun  DC_BETA.BAS 


420  REM  Program  DC_BETA.BAS,  Brian  Kimball,  3/3/92 
500  REM 

510  REM  This  program  calculates  thermistor  constants  from 
520  REM  temperature/voltage  data. 

740  REM 

760  REM  DEFINE  VARIABLES 

780  REM  R  -  voltage  divider  resistor  resistance 

800  REM  TL  -  low  temperature  (thermometer  room  temperature) 

820  REM  VSL  -  source  voltage  at  low  temperature 
840  REM  VTL  -  voltage  across  resistor  R  at  low  temperature 
860  REM  VSLP  “  source  voltage  at  low  temperature  and  low  power 

880  REM  VTLP  -  voltage  across  resistor  R  at  low  temperature  and 

low  power 

900  REM  RTL  -  resistance  across  thermistor  at  low  temperature 
920  REM  ETL  -  voltage  across  thermistor  at  low  temperature 
940  REM  RTLP  -  resistance  across  thermistor  at  low  temperature 
and  low  power 

960  REM  ETLP  -  voltage  across  thermistor  at  low  temperature  and 
low  power 

980  REM  TM  -  mid  temperature  (thermometer  ~36  C) 

1000  REM  VSM  -  source  voltage  at  mid  temperature 

1020  REM  VTM  -  voltage  across  resistor  R  at  mid  temperature 

1040  REM  TAM  -  thermometer  stem  temperature  at  mid  temperature 

1060  REM  RTM  -  thermistor  resistance  at  mid  temperature 

1080  REM  ETM  -  thermistor  voltage  at  mid  temperature 

1100  REM  TH  -  high  temperature  (thermometer  ~50  C) 

1120  REM  VSH  -  source  voltage  at  high  temperature 
1140  REM  VTH  -  voltage  across  resistor  R  at  high  temperature 
1160  REM  TAH  -  thermometer  stem  temperature  at  high  temperature 
1180  REM  RTH  -  thermistor  resistance  at  high  temperature 
1200  REM  ETH  -  thermistor  voltage  at  high  temperature 
1220  REM  PL  -  power  dissipated  in  thermistor  at  low  temperature 
1240  REM  PLP  -  power  dissipated  in  thermistor  at  low  temperature 
low  power 

1260  REM  PM  -  power  dissipated  in  thermistor  at  mid  temperature 

1280  REM  PH  -  power  dissipated  in  thermistor  at  high  temperature 

1300  REM  DCI  -  assumed  value  of  dissipation  constant  (temporary 
variable) 

1320  REM  DC  -  dissipation  constant 

1340  REM  WRT  -  uncertainty  in  thermistor  resistance 
1360  REM  WRTR  -  uncertainty  in  reference  temperature  thermistor 
resistance 

1380  REM  WT  -  uncertainty  in  temperature 

1400  REM  WTr  -  uncertainty  in  reference  temperature 

1420  REM  T  -  temperature  (temporary  variable) 

1440  REM  Tr  -  reference  temperature  (temporary  variable) 

1460  REM  RT  -  thermistor  resistance  (temporary  variable) 
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1480  REM  RTr  -  thermistor  reference  temperature  resistance 
(temporary  variable) 

1500  REM  BLM  -  beta  between  low  and  mid  temperatures 
1520  REM  PI,  P2,  P3,  P4  -  partials  for  calculating  uncertainty 
(temporary  variables) 

1540  REM  WBLM  -  uncertainty  in  BLM 

1560  REM  BMH  -  beta  between  mid  and  high  temperatures 
1580  REM  WBMH  -  uncertainty  in  BMH 

1600  REM  BLH  -  beta  between  low  and  high  temperatures 

1620  REM  WBLH  -  uncertainty  in  BLH 

1640  REM  BAVE  -  average  of  betas 

1660  REM  EXTL  -  experimental  low  temperature 

1680  REM  EXTLP  -  experimental  temperature  at  low  power 

1700  REM  ERRL  -  absolute  error  in  experimental  low  temperature 

1720  REM  EXTH  -  experimental  high  temperature 

1740  REM  ERRH  -  absolute  error  in  experimental  high  temperature 

1760  REM  TLC  -  corrected  actual  thermistor  low  temperature 

1780  REM  TMC  -  corrected  actual  thermistor  mid  temperature 

1800  REM  THC  -  corrected  actual  thermistor  high  temperature 

1820  REM  RTref  -  final  reference  temperature  resistance  at  36 
degrees  C 
1840  REM 
1860  R  =  9998.8 

1862  PRINT  "Certain  variables  within  a  calibration  run  do  not 
change;  it" 

1864  PRINT  "is  only  necessary  to  enter  them  once  per  batch.  If 
you  are  " 

1866  PRINT  "beginning  a  new  batch  answer  'Y'  at  the  prompt. 
Enter  'N'  if" 

1868  PRINT  "you  have  already  defined  'batch'  variables." 

1870  g$  =  "n" 

1872  PRINT 

1874  PRINT  "Are  you  beginning  a  new  calibration  batch?" 

1875  DO  UNTIL  g$  =  "Y"  OR  g$  =  "N" 

1876  INPUT  "Enter  Y/N;  ",  g$ 

1877  LOOP 

1878  IF  g$  =  "Y"  THEN 

1879  PRINT  "YES" 

1880  GOTO  1885 

1881  ELSE 

1882  PRINT  "NO" 

1883  GOTO  2080 

1884  END  IF 

1885  PRINT  "Input  'batch'  variables" 

1890  INPUT  "  TL  =  ",  TL 

1900  INPUT  "VSL  =  ",  VSL 
1920  INPUT  "LOW  POWER  VSL  =  ",  VSLP 
1940  INPUT  "  TM  =  ",  TM 
1960  INPUT  "VSM  =  ",  VSM 
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2000  INPUT  "STEM  TEMP  (M)  =  ",  TAM 
2020  INPUT  "  TH  =  ",  TH 
2040  INPUT  "VSH  =  ",  VSH 
2060  INPUT  "STEM  TEMP  (H)  =  ",  TAH 

2080  INPUT  "VTL  =  ",  VTL 

2100  INPUT  "LOW  POWER  VTL=",  VTLP 

2120  RTL  =  R  *  ((VSL  /  VTL)  -  1) 

2140  ETL  =  (VSL  *  RTL)  /  (RTL  +  R) 

2160  RTLP  =  R  *  ((VSLP  /  VTLP)  -  1) 

2180  ETLP  =  (VSLP  *  RTLP)  /  (RTLP  +  R) 

2200  REM 

2210  INPUT  "VTM  =  ",  VTM 

2220  RTM  =  R  *  ( (VSM  /  VTM)  -  1) 

2240  ETM  =  (VSM  *  RTM)  /  (RTM  +  R) 

2260  REM 

2280  INPUT  "VTH  =  ",  VTH 

2300  RTH  =  R  *  ((VSH  /  VTH)  -  1) 

2320  ETH  =  (VSH  *  RTH)  /  (RTH  +  R) 

2340  REM 

2360  REM  CORRECTION  FOR  EMERGENT  STEM 
2380  TMC  =  TM  +  (.00016  *  111  *  (TM  -  TAM)) 

2400  THC  =  TH  +  (.00016  *  111  *  (TH  -  TAH)) 

2420  REM  CORRECTION  FOR  THERMISTOR  SELF -HEATING 
2440  REM  CALCULATE  POWER  DISSIPATED 
2460  PL  =  ((ETL  ^  2)  /  RTL) 

2480  PLP  =  ((ETLP  "  2)  /  RTLP) 

2500  PM  =  ((ETM  ^  2)  f  RTM) 

2520  PH  =  ((ETH  ^  2)  /  RTH) 

2540  TLC  =  TL  +  273.15 
2560  TMC  =  TMC  +  273.15 
2580  THC  =  THC  +  273.15 
2600  REM  DISSIPATION  CONSTANT 
2620  DCI  =  .0003 
2640  DC  =  .00033 

2650  PRINT  "Iterating  on  dissipation  constant... 

2660  IF  ABS(DC  -  DCI)  <  (.000001)  THEN 

2680  GOTO  4120 

2700  ELSE 

2720  GOTO  2760 

2740  END  IF 

2760  DC  =  (DC  +  DCI)  /  2 

2780  REM  ADD  SELF  HEATING  TERM 

2800  TLC  =  TLC  +  (PL  /  DC) 

2820  PRINT  DC 

2840  TMC  =  TMC  +  (PM  /  DC) 

2860  THC  =  THC  +  (PH  /  DC) 

2880  REM 

2900  WRT  =  3 

2920  WRTr  =  3 

2940  WT  =  .05 
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WTr  =  .05 
REM 

REM  CALCULATE  BETA-LM  AND  BETA  UNCERTAINTY 
T  =  TLC 
Tr  =  TMC 
RT  =  RTL 
RTr  =  RTM 

BLM  =  (T  *  Tr  /  (Tr  -  T))  *  (LOG(RT  /  RTr)) 

REM 

REM  CALCULATE  PARTIALS  FOR  BETA-LM 
PI  =  T  *  Tr  /  (RT  *  (Tr  -  T) ) 

P2  =  T  *  Tr  /  (RTr  *  (Tr  -  T) ) 

P3  =  ((Tr  "  2)  /  ((Tr  -  T)  ^2))  *  LOG(RT  /  RTr) 

P4  =  ((T  ^  2)  /  ((Tr  -  T)  "2))  *  LOG(RT  /  RTr) 

WBLM  =  ((PI  *  WRT)  -  2  +  (P2  *  WRTr)  ^  2  +  (P3  *  WT)  ^  2 
+  (P4  *  WTr)  2)  ^  .5 
REM 

REM  CALCULATE  BETA-MH  AND  BETA  UNCERTAINTY 
T  =  THC 
Tr  =  TMC 
RT  =  RTH 
RTr  =  RTM 

BMH  =  (T  *  Tr  /  (Tr  -  T))  *  (LOG(RT  /  RTr)) 

REM 

REM  CALCULATE  PARTIALS  FOR  BETA-MH 
PI  =  T  *  Tr  /  (RT  *  (Tr  -  T) ) 

P2  =  T  *  Tr  /  (RTr  *  (Tr  -  T) ) 

P3  =  ((Tr  "  2)  /  ((Tr  -  T)  ^2))  *  LOG(RT  /  RTr) 

P4  =  ((T  ^  2)  /  ((Tr  -  T)  "2))  *  LOG(RT  /  RTr) 

WBMH  =  ((PI  *  WRT)  -  2  +  (P2  *  WRTr)  ^  2  +  (P3  *  WT)  ^  2 

+  (P4  *  WTr)  ^2)  ^  .5 

REM 

REM  CALCULATE  BETA-LH  AND  BETA  UNCERTAINTY 
T  =  TLC 
Tr  =  THC 
RT  =  RTL 
RTr  =  RTH 

BLH  =  (T  *  Tr  /  (Tr  -  T) )  *  (LOG(RT  /  RTr)) 

REM 

REM  CALCULATE  PARTIALS  FOR  BETA-LH 
PI  =  T  *  Tr  /  (RT  *  (Tr  -  T) ) 

P2  =  T  *  Tr  /  (RTr  *  (Tr  -  T) ) 

P3  =  ((Tr  "  2)  /  ((Tr  -  T)  "2))  *  LOG(RT  /  RTr) 

P4  =  ((T  "  2)  /  ((Tr  -  T)  "2))  *  LOG(RT  /  RTr) 

REM 

REM  Calculate  temp  at  extremes  using  linear 
approximation. 

REM  Also,  calculate  absolute  error  =  T(actual) -T(calc) 
WBLH  =  ((PI  *  WRT)  ^  2  +  (P2  *  WRTr)  ^  2  +  (P3  *  WT)  "2 
+  (P4  *  WTr)  ^2)  .5 

BAVE  =  (BLM  +  BMH  +  BLH)  /  3 
EXTL  =  (BLH  *  TMC  /  (TMC  *  LOG  (RTL  /  RTM)  +  BLH))  - 
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4000  EXTLP  =  (BLH  *  TMC  /  (TMC  *  LOG(RTLP  /  RTM)  +  BLH) )  - 

273.15 

4020  DCI  =  (PL  -  PLP)  /  (EXTL  -  EXTLP) 

4040  TLC  =  TLC  -  (PL  /  DC) 

4060  TMC  =  TMC  -  (PM  /  DC) 

4080  THC  =  THC  -  (PH  /  DC) 

4100  GOTO  2660 
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DC  =  (DC  +  DCI)  /  2 

TLC  =  TLC  +  (PL  /  DC) 

TMC  =  TMC  +  (PM  /  DC) 

THC  =  THC  +  (PH  /  DC) 

ERRL  =  TLC  -  EXTL  -  273.15 

EXTH  =  (BLH  *  TMC  /  (TMC  *  LOG(RTH  /  RTM)  +  BLH))  -  273.15 

ERRH  =  THC  -  EXTH  -  273.15 

TLCC  =  TLC  -  273.15 

TMCC  =  TMC  -  273.15 

THCC  =  THC  -  273.15 

REM  .  ^  4-0^ 

REM  Calculate  reference  temperature  resistance  at  36 

degrees , Celsius .  ^ » x  v 

RTref  =  RTM  /  (EXP (BLH  *  ((1  /  TMC)  -  (1  /  309.15)))) 

Z1  =  PL  /  DC 
Z2  =  PM  /  DC 
Z3  =  PH  /  DC 
REM 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
g$  =  ''n" 

PRINT  "Do  you  wish  to  continue?” 

DO  UNTIL  g$  =  "Y"  OR  g$  =  ”N” 

INPUT  "Enter  Y/N:  ",  g$ 

LOOP 


•••',  "BETA,  Kelvin”,  "UNCERTAINTY”,  ”RT,  ohms" 

"BIM”,  BLM,  WBLM,  RTL 

"BMH”,  BMH,  WBMH,  RTM 

"BLH”,  BLH,  WBLH,  RTH 

"BAVE",  BAVE 

,  •’TEMP”,  "EXP  TEMP”,  "ERROR” 

"LOW”,  TLCC,  EXTL,  ERRL 
••MED”,  TMCC,  ””,  ”” 

"HIGH”,  THCC,  EXTH,  ERRH 

”RT36  =  ”,  RTref 
"DC  =  ”,  DC 
”SHL”,  "SHM”,  "SHH” 

Zl,  Z2,  Z3 


4827  IF  g$  =  ”Y”  THEN 

4828  PRINT  "YES” 
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4829  goto  1870 

4830  ELSE 

4831  PRINT  "NO” 

4832  GOTO  6000 

4833  END  IF 
6000  END 
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